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Abstract 

Extending our earlier study of nonlinear Bogolyubov-Valatin transforma- 
tions (canonical transformations for fermions) for one fermionic mode, in the 
present paper we perform a thorough study of general (nonlinear) canoni- 
cal transformations for two fermionic modes. We find that the Bogolyubov- 
Valatin group for n = 2 fermionic modes which can be implemented by means 
of unitary SU{2"' = 4) transformations is isomorphic to 50(6;]R)/Z2. The 
investigation touches on a number of subjects. As a novelty from a math- 
ematical point of view, we study the structure of nonlinear basis transfor- 
mations in a Clifford algebra [specifically, in the Clifford algebra C(0,4)] 
entailing (supersymmetric) transformations among multivectors of different 
grades. A prominent algebraic role in this context is being played by bipa- 
ravectors (linear combinations of products of Dirac matrices, quadriquater- 
nions, sedenions) and spin bivectors (antisymmetric complex matrices). The 
studied biparavectors are equivalent to Eddington's ii^-numbers and can be 
understood in terms of the tensor product of two commuting copies of the 
division algebra of quaternions M. From a physical point of view, we present 
a method to diagonalize any arbitrary two-fermion Hamiltonian. Relying on 
Jordan- Wigner transformations for two-spin-^ and single-spin- 1 systems, we 
also study nonlinear spin transformations and the related problem of diago- 
nalizing arbitrary two-spin-^ and single-spin-| Hamiltonians. Finally, from a 
calculational point of view, we pay due attention to explicit parametrizations 
of 317(4:) and 50(6; M) matrices (of respective sizes 4x4 and 6x6) and their 
mutual relation. 
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1 Introduction 



Like canonical transformations in classical mechanics, unitary transformations of 
quantum dynamical degrees of freedom often simplify the dynamical equations, or 
allow to introduce sensible approximation schemes. Such methods have wide-ranging 
applications, from the study of simple systems to many-body problems in solid-state 
or nuclear physics and quantum chemistry, up to the infinite-dimensional systems of 
quantum field theory [1-4]. Linear (unitary) canonical transformations (i.e., trans- 
formations preserving the canonical anticommutation relations (CAR)) for fermions 
have been introduced by Bogolyubov and Valatin (for two fermionic modes) in con- 
nection with the study of the mechanism of superconductivity [5, 10-13]. These 
(linear) Bogolyubov- Valatin transformations have been extended, initially by Bo- 
golyubov and his collaborators [14, 15], [16], Appendix II, p. 123 ([17], p. 116, [18], p. 
679), to involve n fermionic modes [so-called generalized linear Bogolyubov- Valatin 
transformations, see, e.g., [19], "qacTb III [chasf III]/[part III], p. 247 ([20], p. 341, 
[21], p. 259, [22], p. 235)]. Such linear canonical transformations are important from 
a physical as well as from a mathematical point of view. Mathematically, they allow 
to relate quite arbitrary Hamiltonians quadratic in the fermion creation and anni- 
hilation operators to collections of Fermi oscillators whose mathematics is very well 
understood. From a physical point of view, canonical transformations implement 
the concept of quasiparticles in terms of which the physical processes taking place 
can be described and understood in an effective and transparent manner. To apply 
the powerful tool of canonical transformations to the physically interesting class 
of non-quadratic Hamiltonians, however, requires to go beyond linear Bogolyubov- 
Valatin transformations. Certain aspects of nonlinear Bogolyubov- Valatin transfor- 
mations have received some attention over time [23-27], [28], Sec. 2-6, p. 52, [29-54] 
(We disregard here work done within the framework of the coupled-cluster method 
(CCM) [4] which is nonunitary.).^ However, a systematic analytic study of general 
(nonlinear) Bogolyubov- Valatin transformations had not been undertaken until the 
publication of our article [49]. 

In [49] we have initiated a systematic analytic study of general (nonlinear) Bo- 
golyubov- Valatin transformations by studying in detail the prototypical case of one 
fermion mode. The mathematics of one fermion mode is closely relelated to the 
mathematics of the division algebra of quaternions H. In the present paper, we 
now investigate the most general (nonlinear) Bogolyubov- Valatin transformations 
for two fermion modes. The original papers of Bogolyubov [5, 10] and Valatin [11] 
study the application of quasifermion operators constructed linearly (with certain 
real coefficients u and v) from two fermionic modes in the most simple way. Our 
study can, therefore, be considered as the widest possible (take this with a grain of 
salt) generalization of the work of Bogolyubov and Valatin for two fermion modes. 
For example, we will encounter the generalization(s) of the condition ([5], p. 59, (p. 

^Up to the referencing, this paragraph is identical to the introductory paragraph of our earher 
article [49] on the one-mode case of nonlinear Bogolyubov- Valatin transformations. 
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42 of the English transL), [10], p. 796) 

+ = I (1) 

which is the well-known condition for the original (linear) Bogolyubov-Valatin trans- 
formation to be canonical (We display this condition here symbolically without pre- 
cisely defining its meaning.). 

The study of nonlinear Bogolyubov-Valatin transformations requires the study 
of a new mathematical problem not considered previously - nonlinear basis trans- 
formations in a Clifford algebra. This involves transformations among multivectors 
of different grades which can be considered as a mathematical realization of the 
physical concept of a supersymmetry. Furthermore, nonlinear Bogolyubov-Valatin 
transformations of two fermion modes have connections to various concrete phys- 
ical and mathematical problems. Via Jordan- Wigner transformations two-fermion 
systems are related to two-spin- 1 and single-spin- 1 systems which have been stud- 
ied from various points of view in the past and have received broad attention in 
the context of quantum information and computation theory recently. Nonlinear 
Bogolyubov-Valatin transformations correspond to nonlinear spin transformations 
in these systems. As two-fermion systems are closely related to the Clifford algebra 
C(0, 4) (and their cousins of different signature), the study of nonlinear Bogolyubov- 
Valatin transformations is closely connected to Dirac matrices and the Dirac equa- 
tion in general. Consequently, early studies of Dirac matrices and the Dirac equation 
performed by Eddington feature prominently in the list of references accompanying 
the article. And finally, via the relation of two-fermion systems to 4 x 4 matrices the 
study of nonlinear Bogolyubov-Valatin transformations has close connections to the 
representation and parametrization of the group S't/(4), and its related orthogonal 
group SO{Q; M) to which it is the double cover. 

As in our earlier paper [49], the main focus lies on structural and methodical 
aspects of the problem of nonlinear Bogolyubov-Valatin transformations. We start 
in Sec. 2 with a discussion of certain algebraic aspects of a (fixed) set of fermion 
creation and annihilation operators for two fermion modes. In particular, we make 
contact with Clifford algebras (Subsec. 2.1) and the concept of paravectors (Subsec. 
2.2). In the main section 3 of our article we first relate the concept of nonlinear 
Bogolyubov-Valatin transformations to the mathematical concept of nonlinear basis 
transformations in a Clifford algebra (Subsec. 3.1) and then study those in mathe- 
matical detail for the Clifford algebra C(0,4) (Subsec. 3.2). We find that nonlinear 
Bogolyubov-Valatin transformations have the structure of the group 5*0(6; ]R)/Z2 
[eq. (127)]. In section 4 we then apply the obtained insight to the problem of diag- 
onalizing fermion (Subsec. 4.1) and spin (Subsec. 4.2) Hamiltonians. As the whole 
discussion in the sections 3 and 4 is closely linked to the intricate relationship be- 
tween the groups 5*0(6; M) and 5t/(4) in section 5 we study explicitly and in detail 
the relationship between their corresponding 6x6 and 4x4 representation matrices. 
Discussions of various aspects of the present work and their relation with the work 
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of other authors are collected in Sec. 6. Short conclusions are presented in Sec. 7. 
To facilitate the readability of the article, a number of technical details are deferred 
to six Appendices. The article is accompanied by a comprehensive list of references 
which may serve as a guide to the relevant literature. 



2 Algebraic preliminaries 

2.1 Fermion creation and annihilation operators and the as- 
sociated CUfford algebra 

We consider two pairs of fermion creation and annihilation operators d^, d^ {k = 
1,2). We regard the creation operators d^ as the hermitian conjugates of the anni- 
hilation operators d^: d^ = d^ (we will use the latter notation throughout). They 
are acting on the vacuum state |0) the following way 

dfciO) = 0, (2) 

4|0) = \l)k. (3) 

The four- dimensional (complex) Hilbert space C4 (Sometimes, we will refer to it 
as spin space.) of the two- fermion system is spanned by the vacuum state |0), the 
two one-particle states and the two-particle state |2)(2,i) [be aware of the index 
convention in eq. (5)] 

d^dllO) = |2)(2,i) , = -|2)(i,2) , (4) 

(0|did2 = (2|(2,l) = -(2|(i,2) . (5) 

The fermion creation and annihilation operators obey the canonical anticommuta- 
tion relations (CAR), [k, I = 1,2; I4 is the 4x4 unit matrix)^, 

{al,di} = alai + d^ dj = 5ki 1^, (6) 
{dfc,dz} = hk ai + hk ai = 0, (7) 

the latter equation entailing 

{«L«!} = 4 «/ + 4 = • (8) 

It is now useful to consider the following pairs of anti-hermitian operators. 

4'^ = -4'^ =^{ak + 4) (9) 

4^ = -4^^ =ak-al (10) 



^As we will emphasize in this paper on many occasions the concrete aspects of the problem 
under consideration we have decided to indicate the presence of unit operators explicitly. 
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As a consequence of the CAR, eqs. (6)- (8), these operators obey the equation (p, q = 
1,2) 



KUri = -2 5pg5kiU. (11) 
We introduce now the following useful notation: 

C2fc-i = 4^^, k = 1, . . . ,n {= 2) , 

C2k = 4'^ k = l,...,n (=2) . (12) 

The anti-Hermitian operators Ck = —c\ obey the Clifford algebra relation [/c, / = 
l,...,2n (= 4)] 

{cfc,Q} = 2 Ck- ci = 2 Qki = -2 5h I4 , (13) 

gki = -hi , (14) 

where Ck ■ ci denotes the inner product in the Clifford algebra. The operators gen- 
erate the (real) Clifford algebra C(0,4) which is isomorphic to the twofold tensor 
product of the algebra of quaternions H (cf., e.g., [55], Chap. 15, p. 123, [56], Chap. 
16, p. 205). To simplify the further calculations and to minimize sign errors, in the 
following we will raise and lower indices by means of 6ki and not by means of Qki- 



2.2 Some useful algebraic structures in the Clifford algebra 

C(0,4) 

Having related the two-mode pairs of fermion creation and annihilation operators 
to the Clifford algebra C(0,4) we will focus now our attention onto this algebraic 
structure. In this subsection, we will discuss certain algebraic objects in the Clifford 
algebra C(0,4) that we will encounter in the course of the further investigation. 
Their significance will become clear in the further sections only. The somewhat un- 
conventional index convention (including an index -1) to be used in the following has 
been chosen with an eye to possible future generalizations of it to a larger number 
of fermionic modes. 



Paravectors 

Define now 

Ck = Ck , /c = 1, . . . ,4 , (15) 

C(_i) = 2CiC2C3C4 . (16) 

These operators obey the equation {k,l = — 1, 1, . . . , 4) 

Ck,ci \ = -2 6ki l4- (17) 



They are generators of the (non-universal) Chfford algebra C(0,5). Going one step 
further we define 



Cfc = Cfc , /c = -1, 1, . . . ,4 , (18) 

£o = -I4 . (19) 

These operators Ck span the paravector space"' Vq associated with the Clifford alge- 
bra C(0, 5) the same way as the operators span the vector space that is associated 
with the Clifford algebra C(0,4). It should be mentioned here that the choice in 

sign for cq = — 14 is a matter of convenience. Our choice has been suggested by eq. 
(105) further below'*. 

As a consequence of eq. (17) and the ant i- Hermit icity of the operators Ck, the 

operators Ck [i.e., the basis elements of the paravector space associated with the 
Clifford algebra C(0, 5)] obey the relation {k, I = —1, . . . , 4) 

Ci^Ci + CiCk = Ck Ci + Ci C,, = 2 4; U4. (20) 

Note, that eq. (20) has the form of the unitary Hurwitz-Radon matrix problem ([65], 
Subsec. 1.4, p. 25, also see [66], part II, p. 67 and [67]). One can convince oneself 
that the relation 

J i 2 -J 2 -J 2 

r = f "■PI'''' r r r r r (21) 

holds {emnpqrs IS the 6-dimensional totally antisymmetric tensor, e(_i)oi234 = 1; we 
have raised here the indices by means of 6ki and not using gki), and consequently 
also [n^p^q^r^s) 

2 it J ;t ; , ^ 2 

C n Cp C q C s = —i ^npqrs ^ rn ■ (22) 

Biparavectors 

We define now the antisymmetric objects 

■A 2 2^ 2 \ Xt 



Cmim2 Crn2mi ^ \^mi^m2 ^m2^mij ^mim2 ' (23) 

the decomposable/simple biparavectors of the paravector space Vq associated with 
the Clifford algebra C(0,5) (cf. [58], Sec. 3.4, p. 12). They span the 15-dimensional 



■^For paravector space see, e.g., [57], Chap. 13, pp. 254-259, [55], Chap. 16, p. 140, [56], 1. and 
2. ext. ed.. Sec. 19.3, p. 247, [58]. The term paravector has been introduced in [59], p. 22. In a 
context related to our study, paravectors have been used carUer in [60], p. 14, eq. (5.3), (in some 
disguise) in [61], p. 340, eq. (1.2), and in [62], Sec. 2, p. 92; also note [63]. For an early, related 
discussion see [64], in particular, p. 3, eq. (5). 

''A different choice (cq = IL4) would have entailed to modify eq. (53) accordingly. 
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biparavector space A^(^)- To simplify the notation we will often write cm instead of 
[M = (mi, 7712); if capital letters, e.g., M, are used for summations, mi < m2 



is understood]. 

It seems worth noting here that the set of biparavectors cm (supplemented by 
the unit operator I4) is equivalent to the i?- numbers introduced by Eddington [68]^. 
It has been recognized from early on^' that the .E-numbers of Eddington form an 
algebraic system^ which can be obtained by taking the tensor product of two (com- 
muting) copies, Ii, Ji, Ki (IiJi = Ki, etc.) and I2, J2, K2, (I2J2 = K2, etc.), of 
the system of quaternions H (cf. also the comment at the end of Subsec. 2.1). To 
be more specific, one possible (symmetric) choice for the generators of the Clifford 
algebra C(0,4) is (for a different choice see Appendix D) 



Cl 



C2 



C3 



Cl 



C2 



C3 



C 4 — C4 



1 

71 

1 

1 
1 

7! 



[I4 - i K1K2] ii 

[l4-Z K1K2] Jl 

[14 + 1 K1K2] I2 
[14 + 1 K1K2] J2 



1 

71 

1 

71 

1 

7! 
1 

71 



[Ii - ^ J1K2] , 

[Jl + ^ I1K2] , 

[I2 + i K1J2] , 

[J2 - ^ K1I2] , 



and, consequently. 



•(-1) 



t K1K2 



(24) 
(25) 
(26) 
(27) 

(28) 



Inversely, the two commuting copies of the system of quaternions can be given in 
terms of the generators of the Clifford algebra C(0,4) by means of the following 
equations. 



Ii = ^ 



1 

V2 



Jl = ^ 



1 

V2 



I4 + C(_i) 



I4 + C(„i) 



Cl 



C2 



1 

7! 

1 

7! 



Coi + C(_i)i 



C02 



■(-1)2 



— = [ci + 7 C2C3C4] (29) 
v2 



-j=[c2-i C1C3C4] (30) 



V2 



^Also see [69-71]; a systematic presentation is given in [72], Part I, Chaps. 2-4, pp. 20-61, 
[73], Chaps. VI/VII, pp. 106-158, and a discussion of Eddington's iJ-numbcrs from a modern 
perspective can be found in [74]. 

6 [75-79], [72], §3.8, p. 47, [80], [73], §60, p. 120, [81, 82]^. 

^One can convince oneself that the biparavector space /\ (Vg) supplemented by the unit operator 
IL4 spans the operator space of the Clifford algebra C(0,4). Initially, biparavectors (supplemented 
by the unit operator) have been named ^ sedenions' ([75], p. 105, footnote 3) or ^ quadriquaternions' 
[83]. Note, that these sedenions are different from the weird algebraic objects presently called 
'sedenions' and obtained from the octonions in the following step of the Cayley-Dickson process. 
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Ki 


= C1C2 




- C12 


= C1C2 












h 


1 

V ^ 




-C(_i) 


C3 = 


1 

V 


Cos 


-C(_i)3 


1 

71 

V ^ 


[C3 - « C1C2C4 


J2 


1 


I4 


-C(-i) 


C4 = 


1 

"71 


C04 


- C(_i)4 


1 

" 71 


[C4 + 2 C1C2C3 


K2 


= C3C4 




~ C34 


= C3C4 













(31) 



(34) 

As the system of quaternions I, J, K is related to the group SU(2) and, therefore, 
to spin variables S^, S'^, {S^ = ^I/2, = iJ/2, S'^ = iK/2), a quadriquater- 

nion representation of the biparavectors cm is related to a two-spin-| system'^. For 
example, this fact is relied on in the product operator formalism widely used in the 
field of nuclear magnetic resonance (NMR) [85-87]. Two-spin- 1 systems will be dis- 
cussed in Subsection 4.2.1. Another approach relating quaternions to Eddington's 
^'-numbers has been put forward by Conway [88] (also see [89], Sec. 9, p. 48). 



Algebra of biparavectors 

Using the notation introduced in eq. (23), the eqs. (21), (22) can be transformed to 
read^ 

Cm = ^ cp cq (35) 



and 

Cp Cq = — 6pQ I4 — 5p^q^ C(_p2q,^) — 5p^q^ C(^p^q^) 

- i epQ^ Cn . (36) 

In the last equation we have partially lifted the condition Pi ^ P2 ^ Qi ^ (I2, it still 
applies Pi 7^ P2, Qi Q2- Equation (36) can be found in a number of studies^°. The 

®In fact, the equations (24)-(34) are a modified version of a Jordan- Wigner transformation used 
in [84], p. 13835, above of eq. (5). 

^Related equations can be found in: 

[90], p. 203, eq. (2), [91] p. 756, eq. (8), [92], p. 27, eq. (14), 

[94], p. 359, eq. (2.8), [95], p. 118, eq. (5), and [96], p. 205, eq. (9). 

[97], p. 656, §8, eqs. (8.6), (8.7) (here, a five-dimensional formulation is still being used), 
[98], p. 32 (reprint # 1: p. 50, reprint # 2: p. 943), eq. (9), 

[90], p. 203, eq. (2), [91], p. 756, eq. (7), [92], p. 27, eq. (12), 

[99], p. 3121, eq. (F), [100], p. 136, eqs. (2.4), (2.10) [101], p. 127, eqs. (2), (3), 
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commutator [cp,Cq\ reads^^: 

[cp,Cq] = CpCQ - CqCp = - 2 C (pjqj) - 2 (5p252 C(piqi) 

while the anticommutator {c p, cq} is given by^^: 

{cp,Cq} = CpCQ + Cq Cp 

= {- "^SpQ +2 6p,q, 6p,q,) I4 - 2i epQ^ cn ■ (38) 

Equation (37) displays the structure of the Lie algebra of su(4) ~ so(6). Note 
that this Lie algebra (i.e., its generators) is realized here in terms of products of 
one to four fermion creation and annihilation operators. This representation of the 
Lie algebra su(4) ^ so(6) has been discussed earlier in [112, 113] (In general, for 
n fermionic modes the Lie algebra generated by all possible products of fermion 
creation and annihilation operators is su(2") [114], Sec. 5.2, p. 21, [115].). Most 
discussions, in the literature, of Lie algebras in terms of fermion creation and an- 
nihilation operators rely on their bilinear combinations only (see, e.g., [116], [117], 
Subsecs. 1.59, 1.60, pp. 88-93) or, at best, supplement these combinations by the 
fermion creation and annihilation operators themselves [31, 35, 118], [119], App. 
1, p. 919 (p. 548 of the reprint), [40, 121], [122], p. 907, [49]. Cases that go be- 
yond this linear /bilinear framework have been discussed in [114], Sec. 5.2, p. 21, 
[112, 113, 115, 123, 124]. Finally, we mention that certain graphical representations 
of the commutator/anticommutator relations (37), (38) have appeared in the liter- 
ature [125-130]. To relate the works of Saniga and collaborators [126-129] and Rau 
[130] to the equations (37), (38) one must rely on the representation of the bipara- 

vectors cm in terms of the (twofold) tensor product of quaternions H discussed at 



[94], p. 359, eq. (2.6), [102], p. 183, eq. (8a), [95], p. 118, eq. (4), 

[103], p. 140, eq. (2.15), [104], p. 2244, eq. (V), [96], p. 204, eq. (1), 

[105], Appendix C, p. 105021-19, eq. (C.13). 

^^The commutator can be found in the literature in: 

[106] p. 97, eq. (1"), [91], p. 756, eq. (7), [92], p. 27, eq. (12), 

[107], p. 354, eq. (2), [108], p. B841, eq. (3.2), [lOO], p. 136, eq. (2.11), 

[109], Subsec. 2d, p. 21, [94], p. 360, eq. (3.2), [102], p. 182, eq. (6), 

[110], p. 368 (p. 1212 of the English transl.), below from eq. (7), 

[95], p. 118, eq. (3), [104], p. 2244, eq. (Ill), [96], p. 204, eq. (3), 

[111], Sec. 5.3, p. 168, [105], Appendix C, p. 105021-19, eq. (C.12). 

^^The anticommutator can be found in the literature in: 
[98], p. 33 (reprint # 1: p. 50, reprint # 2: p. 944), eq. (10a), 

[91], p. 756, eq. (7), [92], p. 27, eq. (12), [94], p. 360, eq. (3.3), 

[102], p. 182, eq. (7), [95], p. 122, eq. (40), [61], p. 341, eq. (1.11). 
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the end of Appendix D. 



3 Nonlinear Bogolyubov-Valatin transformations 

3.1 Clifford algebra formulation of canonical fermion trans- 
formations 

3.1.1 The ansatz for nonlinear Bogolyubov-Valatin transformations 

We proceed now to writing down an ansatz for the most general Bogolyubov-Valatin 
transformation for two fermionic modes. In view of the eqs. (6), (7), the new fermion 
annihilation and creation operators bk, hi, read 

bk = B,.({A};{a}) = U a, 

- \(0|0) , \{1|0) f,^ , \{2|0) p^^ , x(0|l) . ,(0|2) . 

, (1,2|0) .fat , p^ip^ , ),(1|2) p^U 

+ -^f ^2^1 + -^f 4^2 + -^1°'^'^^ CLid2 
, (1|1,2) .f. . . (2|1,2) .f;; 

+ A^^'^'^'^^ a{dlaia2 , (39) 
6l = B,({A};{a})t = UalUt 



, (010) ^ , (110) . , (2|0) . , (0|1) .f , (0|2) 

-^fc -"-4 + -^k '^1 + -^fc '^2 + «1 + ^fc (^2 



, (1,2|0) , x(l|l) ;^t;i , x(l|2) p^U 

A/j Ciici2 + A^ a^^fli + A^ a2ai 



+ Af ata2 + Al^l^) 4a2 - aI^I^'^^ 44 



— A[. ' a|aia2 — -^i ' ^ 02^102 



A^ ' ' o.|ci2^i — A^ ' ' ci|ci2^2 



+ aJ^^'^'^'^^ al4«ia2 . (40) 

We assume the coefficients to be complex numbers: A*^"'^ G C. {A} denotes the set 
of all 16 coefficients A^ ' ^ and {a} the set of the creation and annihilation operators. 
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U = U ({A}; {a}) is an unitary operator belonging to the group 5'f/(4) that im- 
plements the nonlinear Bogolyubov-Valatin transformation. The new annihilation 
operators 6^ act on the new vacuum state 

|0') = |0,{A}) = U({A};{a})|0) (41) 

the usual way 

h\0') = 0. (42) 

The new fermion creation and annihilation operators also must obey the canonical 
anticommutation relations [CAR, cf. eqs. (6), (7)], {k,l = 1,2), 

{blk} = bik + h b] = 5ki U , (43) 
{bk^bi} = bkk + bk k = 0, (44) 

the latter equation entailing 

{^,^!} = im + hm = 0. (45) 

Before proceeding further, it is useful to consider the expression for the trace of 
any annihilation (or creation) operator. 

tr bk 

= {0'\h\0') + (l'|i&fc|l')i + (l'|2^fc|l')2 + (2'|(i,2) &fc|2')(i,2) = (46) 
The trace has to vanish. From eq. (39), for eq. (46) it then follows: 

tr h = 4Af °) + 2Ai^l^) + 2A1^I^) - A^^l^'^^ = (47) 

(In the one-mode case studied earlier [49], the corresponding equation is eq. (13), p. 
10247.). 

To facilitate the further discussion, we now introduce the antihermitian operators 

^i^^ = -r =^(;h+bi) 

= Re 4^1°) a? + Re 4°!^) a? + Re 4^1°^ a™ + Re 4°!^) a^^ 
+ Re 4^W a?a? + Re 4^-^l°) a^a^^ + Re 4^1^^ a?a^^^ 
+ Re 4^1^) aS^W + Re 4°!^'^^ a?^a^^ + Re 4^1^^ aWaf 
+ Re 4^'^W .aWaS^W + 4^1^'^^ ^a?df^a^^ 
+ Re 4^'^!^) zaWaWa'^l + Re 4^1^'^^ ta?^a^^a^^ 
+ Re 4^'^l^'^) zaWaflaWa|l (48) 
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and 

= -r =h-hi 

= Im 4^1°) + Im af + Im nf'^ a™ + Im nf^^ af 
+ Im aflafl + Im 4^'^l°) a^aW + Im 4^1^^ a^'^af 

+ Im 4^W ^ ^(011,2) .[2].[2] ^ ^(2|2) .[1] .[2] 

+ Im 4^'^'^'' ^2^' + Im 4^'^'^^ mj^^^df^flg^^ 

+ Im 4^'^'^'' ^tir'a2^'a2^' + Im 4^'^'^^ icL{^a^2^'2^ 
+ Im 4^'^!^'^) zaWaf]aWa|l . (49) 

With some hindsight, the primary ordering of the operators on the r.h.s. is being 
done in accordance with the mode number. In deriving the eqs. (48), (49) we have 
taken into account eq. (47). The explicit relation of the coefficients in the above 
two equations to the coefficients in the eqs. (39), (40) is given in Appendix A. As a 
consequence of the CAR, eqs. (43)- (45), the operators 6^' must obey the equation: 

{^Uf^} = -2(5p,4zl4. (50) 

3.1.2 Clifford algebra counterpart of nonlinear Bogolyubov-Valatin 
transformations 

In analogy to eq. (12), we introduce the following useful notation: 

4fc-i = , k = l,...,n (=2) , 

d2k = , k = l,...,n (=2) . (51) 

The operators dk must obey the C(0,2n = 4) Clifford algebra relation [k,l = 
l,...,2n (= 4)] 

|4,cli} = 2 dk-di = 2 Qki I4 = -2 Ski I4 • (52) 

This relation is the Clifford algebra analogue of the canonical anticommutation rela- 
tions, (43)- (45), for fermions. The study of nonlinear Bogolyubov-Valatin transfor- 
mations can consequently be performed by studying nonlinear basis transformations 
in the Clifford algebra C(0,4). 

To prepare ourselves for the further investigation, the equations (48), (49) now 
can be compactly written as [The factor of i in the last two terms can also be 
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deduced, in principle, from the condition = — see eqs. (48), (49).]: 



; _ [1] (m) ^ 1 [2] (m,n) ^ ^ _|_ J. [3] (rn,n,p) .... 



~^ ^\^k ^ ^ ^ ^mCnCpCq (53) 

= xL'' ^"^^ + ^xL" ^"'"^ A c„ + ixL'^ ^'"'"'^^ ^ c„ A c„ A c, 

+ ^X?^ ^'"'"'"''^ ^ a„ A c„ A c, A a, (54) 

, [2] (1,2) ... [2] (1,3) ... [2] (1,4) . . 
+ Xk C1C2 + XV C1C3 + XV C1C4 

, [2] (2,3) ... [2] (2,4) . . , [2] (3,4) . . 
+ Xk C2C3 + XV C2C4 + Xk C3C4 

, [3] (1,2,3) .... , [3] (1,2,4) .... 
+ Xk ^ C1C2C3 + Xk ^ C1C2C4 

, [3] (1,3,4) .... [3] (2,3,4) .... 

+ Xk ^ C1C3C4 + Xfc ^ C2C3C4 

, [4] (1,2,3,4) ..... 
+ Xk ^ C1C2C3C4 (55) 

The tensors x|,^' (™'")^ ^[3] {"i:".p)^ {m,n,p,q) ^^^j ^j^^ j^^^ three ten- 

sors are totally antisymmetric with respect to their upper indices. For the relation 
of the tensors xL'^ "* to the coefficients k^'"-* see Appendix A. 

In principle, already here one could go one step further and define the operators 
dk in terms of biparavectors (with some coefficients Xk *^)- 

4 = - 4 = - Xk Cm • (56) 

However, this step can much better be understood in course of the further investi- 
gation and we postpone it, therefore, to Part II of Subsec. 3.2.3, eq. (105). 



3.2 Nonlinear basis transformations in the Clifford algebra 

C(0,4) 

3.2.1 Anticommutators — Conditions for canonical transformations 

From the eqs. (53), (54), we find for the anticommutator, eq. (52), 



|4,(i;| = dkdi + didk =^ 2dk-di 
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Ji] .Ji] M 



[2] , [2] (m,n) 



(m,n) 



21 '^fc (m,n) 



3! (™>"->p) 



[3] {m,n,p) 



J4] 



[4] (m,n,p,q) 



— I 



[2] [3] (m,n,g) [3] (m,n,g) [2] 

X-k {m,n) M "r Afc A/ (m,n) 



[1] [3] (m,p,g) [3] (m,p,g) [1] 

Afc (m) A/ "T Afc A/ (m) 



, 1 [2] [4] (m,n,p,g) 1 [4] {m,n,p,q) [2] 

"""o -^fe (m,n) M 9 ^' ' (m,n) 



Cp A Cqr 



[1] (n) [2] (p,q) [2] (p,g) [1] (n) 

/(-k A,/ ~r Afc A/ 



A Cp A Cq 



+ 2 ^fc 



[2] (m.n) ^ [2] (p,q) 



Cyji A c^i, A Cp A Cg . 



(57) 



The above result has been derived by means of the eqs. (B.1)-(B.4) given in Appendix 
B [In the product dkdi, use for the first factor the representation (53) and for the 
second factor (54).]. Now, imposing on the anticommutator (57) the Chfford algebra 
analogue of the canonical anticommutation relation (CAR) 



|4,rfij = 2 Qki I4 = -2 5ki 
the equations (59)- (63) follow: 



(5^ 



scalar part of eq. (58): 



^[1] ^[1] {r,^) £ [2] [2] (m,n) 

Afc (m) Ai I 21 X-k (m,n) M 



J3] 



3! 



Xk (m,n,p) Xl 



[3] {m,n,p) 



^1 Xk {m,n,p,q) Xl 



[4] (m,n,p,q) 



^[1] (1) ^[1] (1) , ^[1] (2) [1] (2) [1] (3) [1] (3) [1] (4) [1] (4) 

Afc Xl ' Xk Xl ~r Xk Xl ~r Afc Xl 

[2] (1,2) [2] (1,2) [2] (1,3) [2] (1,3) [2] (1,4) [2] (1,4) 

~r Afc Ai Afc Xl ~r Afc Ai 

, [2] (2,3) [2] (2,3) , [2] (2,4) [2] (2,4) , [2] (3,4) [2] (3,4) 

+ Xk Xl + Xk Xl + Xk Xl 

[3] (1,2,3) [3] (1,2,3) [3] (1,2,4) [3] (1,2,4) [3] (1,3,4) [3] (1,3,4) 

' Xk Xl ~r Afc Xl ~r Afc Xl 



, [3] (2,3,4) [3] (2,3,4) , [4] (1,2,3,4) [4] (1,2,3,4) 
+ Xk Xl + Xk Xl 



6k. 



kl 



(59) 
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(The above equation is the orthogonahty/orthonormahty condition for 4 vec- 
tors in a 15- dimensional Euchdean vector space.), 

vector part of eq. (58): 



■ [2] [3] im,n,q) [3] {m,n,q) [2] 

■^k {m,n) {m,n) 



(60) 



bivector part of eq. (58): 



[1] [3] (m,p,q) [3] (m,p,<?) [1] 

Afc (m) A/ "T Afc A/ (m) 



.[2] 



~'~o ^fc (m,n) 



[4] {m,n,p,q) 



+ 2 ^fc 



[4] {m.n.p.g) [2] 



/ (m,n) 



Cp A c, = , (61) 



trivector part of eq. (58): 



" [1] (n) [2] (p,g) [2] (p,g) [1] (n) 

Afc a; "t Afc a; 



A Cp A 



(62) 



quadrivector part of eq. (58): 



[2] (m,n) [2] {p,g) ^ a " a " a " 
Xk Xl A C„ A Cp A Cq 







(63) 



These equations have to be studied in order to investigate basis transformations in 
the Chfford algebra C (0, 4) and, consequently, nonlinear Bogolyubov-Valatin trans- 
formations for two fermionic modes. 

For the sake of completeness, here we also give the expression for the commutator 
of the operators dk (which, incidentally, also defines their wedge product). 



dk, di 



dud, - didk = 2 dk A d, 



[2] {m,q) [I] _ [1] [2] (m,g) j_ [4] {m,n,p,q) [3] 

Xk Xl (m) Afc (m) Xl ' gj Afc A; {m,n,p) 



3! 



Xk Xl 



[4] {m,n,p,q) 



+ 2 



^[1] ^[1] (9) , ^[2] [2] {m,g) 

Xk (p) Xl ' Afc (m,p) Xl 



, 1 [3] [3] {m,n,q) 

2 ' 



Ac„ 



1 [4] {m.,n,p,q) [1] 
q Afc Xl (m) 



1 ^[1] ^[4] {m,n,p,q) 
q Afc (m) Xl 



17 



+ 



[2] (r,n) r [3] (s,p,g) _ [3] {s,p,q) . [ 
Afc "rs Xl Ak "rs Xl 



[2] {r,n) 



Cfi A Cp A 



+ 



^[1] M^[3] (n,P,g) 



[3] {n,p,q) [I] (m) 



A c„ A Cp A Cg (64) 



At the end of this subsection, a conceptual comment is in order. The first two 
hues of eq. (57) define the inner product in the new set of vectors dk- As the vectors 
dk are not hnear combinations of the original vectors Ck, the vector space within 
the Clifford algebra C(0,4) spanned by the vectors dk differs from the vector space 
spanned by the vectors Ck ^'^ Consequently, also the wedge products built over the 
respective vector spaces differ conceptually, as do concepts related to them such as 
the grade of a multivector. Consequently, the first line of eq. (64) defines the wedge 

product related to the vector space spanned by the vectors dk (We indicate this fact 

d 

by means of the letter 'd' atop of the wedge symbol: A.). To see this explicitly 
compare eq. (64) with the following equation which has been evaluated using the 
linearity and associativity of the wedge product related to the vectors Ck- 



dk A di 



[1] [1] (<?) c f A c 
Xk i'p) Xl ^ ^ 



+ 



2! 



[1] (n) [2] 
Xk Xl 



m (p,q) [1] {n) 



Cn A Cp A Cq 



+ 



^ , fl] (m)^ [3] (n,p,g) « ^ [3] (n,p,g)^il] (m) 



3! 



T Xk Xl 



3! 



Xk Xl 



-L [2] (m,n) [2] (p,g) 
^ Xk Xl 



CjYi A A Cp A Cq 



(65) 



Contrary to what naively one might expect, the product dk A di fails to be antisym- 
metric with respect to the indices k and / unless the condition dk ■ di = —Ski I4 [eq. 
(58)] is imposed. 



3.2.2 Solution of the Clifford algebra analogue of the CAR — Infinitesi- 
mal method 

We are now going to solve the eqs. (59)-(63) by means of two different meth- 
ods, an infinitesimal and a global one. We start with the infinitesimal method 
in the vicinity of the identical Bogolyubov-Valatin transformation. We apply the 

relations Xk (m) ^ '^^'^ + '^Xk (m)' Xk (m,n) ^ ^^fc (m,n)' Xk {m,n,p) ^ ^^fc (m,n,p)' 

xt\m,n,p,q) = ^xt\m,n,p,q)^ whcre all quantities preceeded by A are infinitesimally 

^■^Certain special types of such nonlinear basis transformations in a Clifford algebras have been 
considered previously in [46, 131, 132] and in generality for the Clifford algebra C(0, 2) in [49]. For 
related aspects see [133]. 
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small. Neglecting higher order terms the eqs. (59), (61), (62) yield the following 
conditions, respectively: 

Axl^o + = ' (66) 

AxflU) + Axi^\;„, = . (68) 

Here, Ax|^^^*^ = eimpq ^Xfe ' ^^''^^ {^impq is the completely antisymmetric tensor in 
four dimensions with €1234 = I-)- To arrive at this form of the eq. (68) we have made 
use of the relation 

^mnpq ^rst ^mr^ns^pt ^ms^nt^pr ~l~ ^mt^nr^ps 

^ms^nr^pt ^mr^nt^ps ^mt^ns^pr • (^9) 

The eqs. (60) and (63) do not give rise to any condition for the infinitesimal param- 
eters. From the eqs. (66)-(68) we can conclude that Ax^^\i^, Ax{f^*;,„), ^xt\m,n,i) 
are objects which are completely antisymmetric in all of their lower indices. Conse- 
quently, we find that the infinitesimal (nonlinear) Bogolyubov-Valatin transforma- 
tions depend on 15 parameters: 6 parameters given by Axj,^'^^^, 4 parameters given 

by Ax^^^*^i^y 1 parameter given by the totally antisymmetric object Axt\mni) 
4 more parameters given by ^X^k\mnpq) (The latter object does not receive any 
restrictions within our infinitesimal consideration.). The number of parameters for 
each object is in agreement with the global solution to be discussed further below 
[cf. eqs. (80), (74), (84), (89)] and the total number of 15 parameters stands in cor- 
respondence to insight coming from the implementation of the Bogolyubov-Valatin 
transformation by means of an unitary transformation U ({A}; c) [cf. eq. (39)] 

4 = B({A};c) = U({A};c) CfcU({A};c)^ (70) 

U ({A}; c) is for two (n = 2) fermionic modes an element of SU(2"' = 4) which is a 
15-parametric group and the double cover of the group 5*0(6; M). 



3.2.3 Solution of the Clifford algebra analogue of the CAR — Global 
method 

I. Structural analysis of the Clifford algebra analogue of the CAR 

We turn our attention now to the global study of the eqs. (59)-(63). The best strat- 
egy appears to be to start with the consideration of the eq. (63), the quadrivector 
part of these equations. The reason for this lies in the fact that in difference to the 
other eqs. (59)- (62), eq. (63) depends on one type of coefficients [x^^ ("^'"■)^ j ^^le 
bivector coefficients in eq. (53)] only. 
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(A) Quadrivector part 

For the case /c = /, eq. (63) reads 



Xfr^xf = (71) 

where we have apphed the notation xt'' — xlf' <^m A c„ for the bivector x^^\ 
Equivalently, this condition can be formulated as ([134], §7, p. 101, eq. (3) [p. 2250, 
eq. (7.3) of the Enghsh transl.]) 

(Xfc' jXL^) = ^'^'^^'^x[\m,n)Xi\p^q) = 0. (72) 

Here, {A, B) denotes the scalar product of the bivectors A, B. For its definition 
see, e.g., [134], § 1, p. 85, eq. (3) [eq. (1.3), p. 2240 of the English translation]. 
Xfc'* denotes the bivector that is related to the Hodge dual of xf^: xf^* '•'"'"^ = 
gmnpg^p]^^ Finally, a third, equivalent formulation for the condition (71) is: 

PfxL'^ = (73) 

where Pf denotes the Pfaffian of the matrix x^^ . The eqs. (71)-(73) each are the 
condition for the decomposability of the bivector xf^ Consequently, we can write 
for the matrix x^^ ^"^'^^ and the corresponding bivector [For the last part of eq. (76) 
use eq. (B.5).] 

^[2] im,n) ^ ^(^) ^(n) _ ^(rn) ^ ^^^^ 

XL'^ = 2f3A'Hk, (75) 

= 2 /3 A = 2 - 2 HkP^ I4 . (76) 

Here, lik = ^^"^ c^. Tik, P denote row vectors (1x4 matrices) The ansatz (74), 
(76) immediately also solves eq. (63) for k 1. 

xL^Axf = (77) 



1^ [135], Vol. I, p. 26 (p. 19 of the English transl.), 

[136], all eds.. Chap. II, § 6, p. 27, 

[137], p. 25, Exercise 5.1, 

[138], Vol. 1, Chap. 3, end of § 5, p. 165, 

[139], p. 69, 

[134], p. 101 (p. 2250 of the English transl.), 

[140], Chap. 8, end of § 4, HpHMep 3 [Primer 3]/[Exaniple 3], p. 341, 2. and 3. ed. (p. 324 of the 
English transl.), 

[141], Vol. 2 (HacTb II: JlnHeiiHafl Ajire6pa [Chasf II: Lineinaya Algebra]), rji. 6 [Chap. 6], § 5, 
pp. 293-295. 

^^Their elements will finally turn out to be real and can, for simplicity, immediately be chosen 
to be real in view of the fact that ' G M. 
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(B) Trivector part 

Eq. (62) can be written as follows. 

Introducing the ansatz (76) into eq. (78), it can be transformed to read 



(3^ 



X^i^A-Hi - -HkAxPl = (79) 
and its solution can immediately be read off to be {'H^i^\ are real numbers) 

= (3io) _ uf /3 . (81) 
The significance of the chosen notation will become evident in a moment. 



(C) Vector part 

We can now also study eq. (60). Introducing the Hodge dual of xj?^ by writing 



^[3] {n^,n,p) ^ ^„,„p, ^[3]*^^ ^^^^^ 

[2] mnpq [3]* , [3]* mnpq [2] — fl ('S9'l 

Afc (m,n) ^ A/ (g) -T Xfc (g) fc A/ (m,n) ~ ^ ' 

Operating on it with the expression erstp c^'Ac^Ac^ and summing over repeated 
indices, eq. (82) can be brought to the form [taking into account eq. (69)] 

A xpl* + xr A = . (83) 

This equation is analogous to eq. (78), and immediately its solution can be found 
to be ini~^\ are real numbers, the sign is chosen with hindsight) 

xt = - {^^-'^ -Hk - nir'^^) , (85) 

or 

xL'^ ^'"'"'"^ = - H, (,) - ni~'^P^,)) , (86) 

xf = - n - (87) 



where we have applied the notation x^ ' = Xfc^ ('"'"•p) Cm A c„ A Cp for the trivector 
xlf , and nl ^""'"'^^ = Uk (<?), /3* (""'"'^^ = e™ /3(g). Again, the significance of 

the chosen notation will become evident in a moment. 
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(D) Bivector part 

We turn our attention now to eq. (61). Introducing the Hodge dual of xlf' by writing 
^mnpq ^W* making use of the eqs. (74), (80), (62) it reads: 



[4] {m,n,p,q) 



- "Hk (n) Xf^* - Xk^* (^(rn) "H/ („) 

And the solution of this equation is found to be: 



[4]* 



[4]* 
Xk 

[4] {m,n,p,q) 
Xk 



mnpq 



(E) Scalar part 

The remaining task consists in studying eq. (59). Inserting the eqs. 
(90) into it, it can be transformed to read: 



= . (88) 

(89) 
(90) 

50), (74), (86), 



n=-l 



(n) ^(n) 



.m=— 1 



4 

E 

n=~l 



^(n) ^(n) 



4« • 



(91) 



This equation obviates the usefulness of the notation chosen earher in the subpara- 
graphs (A)-(D). From now on we denote by the symbol "H the 4x6 matrix T-i with 
row number k and column number m matrix elements 'H^,™'' and by (3 the row vector 
(1x6 matrix) fi. Eq. (91) can then compactly be written as 

- (-H/?^) {m^) = I4 . (92) 

It is immediately obvious from eq. (92) that it is invariant under (S')0(6;]R) trans- 
formations [H' = HL, = L is a 6 X 6 matrix with L e (^)0(6;M)]. 
Consequently, any two-mode nonlinear Bogolyubov-Valatin transformation can be 
parametrized by means of the (S')0(6; M) transformed parameters of any particular 
(fixed) Bogolyubov-Valatin transformation. The most natural starting point is the 
identical Bogolyubov-Valatin transformation = c^, i.e., in view of eq. (80) we can 
choose 



n 



/3 



(m) 



(3i 



(0,1,0,0,0,0) 



(m) 



6. 



k ' 



k, m 



n 



(-1) 



H 



(0) 



. 



(93) 
(94) 

(95) 
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We immediately find that always 

= , (96) 



/3f3^ = = 1 (97) 

applies. Consequently, eq. (92) finally reads 

nn^ = I4 • (98) 

We can now go one step further in streamlining the notation by writing /3 = "Hq 
(i.e., = T^o"*^). From now on we denote by the symbol "H the 5x6 matrix T-L 
with row number k = 0, . . . ,4 and column number m = — 1, . . . , 4 matrix elements 
-^(m)^ The equations (96)-(98) can then compactly be written in a single line as 

= I5 • (99) 

Any nonlinear Bogolyubov-Valatin transformation can conveniently be parametrized 
by means of a {S)0{6; M) matrix L according to the equation 

H = Hi L . (100) 

The above structural analysis of the Clifford algebra analogue of the CAR has 
now paved the way for a systematic discussion of nonlinear basis transformations in 
the Clifford algebra C(0,4) (and the nonlinear Bogolyubov-Valatin transformations 
related to them). 



II. Biparavector structure of basis transformations in the Clifford algebra 

C(0,4) 

To structure the results obtained in the previous subsection, let us define the 
antisymmetric 6x6 matrix xok with matrix elements XoT^'™^^ means of the 
equation[M = (mi, 7712)] 



(mi,m2) 
AOfc 

or, symbolically. 



Xok 



M 



nj{rni) nj(m2) 
^0 



Xok 



[4]* 
Xk 

[3]* 
Xk 



T 



[4]* 
-Xk 



x'^ 



nj{m2) njirni) 
^0 ^k 



C2(W)„,""™ (101) 



T 
Xok 



J 



(102) 



2noAnk = 2 [CHiAni) 02(1)] 



Ok 



Here, Xk\ xt^* denote the row vectors (1x6 matrices) xlf'* with components 
Xfc^ Xfc'* respectively, and C/(L) denotes the order / compound matrix of 
the matrix L (cf. Appendix C). Writing eq. (53) as 



[1] J3]* 



(103) 



dk 



Xk 



[1] (m) 



'^m + 2] Xk 



[2] (m,n) 



Crn.Cn 



[3]* (g) J ^ , [4]* J 



(104) 
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one can convince oneself that the following representation applies [here, M 
(mi, 7722), mi < m2]. 



du 



di 



1 $t I 

2 c Xofc c 



-2 {-HoAnkf' c 



M 



(105) 
(106) 



We use the notation 



Co 

Cl 
C2 
C3 
V £4 / 



/ iCiC2C3C4 \ 
-I4 
Cl 
C2 

C3 

V C4 / 



(107) 



;t xt xt it it ;t 
C_i, Cq, Cl, C2, C3, C4 



^-ZCiC2C3C4, -I4, -Cl, -C2, -C3, -C4) 



(108) 



Cm used in eq. (105) is a biparavector [cf. eq. (23)]. Taking into account the eqs. 
(93)-(95), (103), eq. (106) can also be written as 



dk 



Ok Cm 



(109) 



This equation describes the structure of nonlinear basis transformations in the Clif- 
ford algebra C(0,4). However, as these transformations must have a group theoret- 
ical structure the eqs. (106), (109) fail to reflect this structural, group theoretical, 
aspect because the new (transformed) generators dk of the Clifford algebra C(0,4) 

are given in terms of other (however, related) objects - the biparavectors cm- 

We will now derive a more symmetric representation of the nonlinear basis trans- 
formations in the Clifford algebra C(0, 4). Taking into account eq. (23), we find (with 
some hindsight) for eq. (106) 



di 



d 



Ok 



c?Q dk — dj^ do 



2 {noA-Hk 

t 2 

— Hn 'Hk 



Cm 

it 

l-Lk ^0 



M 



Ok 



Cm 



(110) 
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Here, Hk = 'HkC = Z]m=-i ^fc"^ is a paravector , = —21-6^^ I4 — iik- From 
the eqs. (20), (99) we recognize that for the paravectors "H^ apphes (/c, / = 0, . . . , 4) 

J t ; X t X X J t X X t 

T-Lf.'Hi + TiiTik = 'Hk'Hi + T-iiHj^ = 2 (5^^ I4 • (m) 
Consequently, we obtain 

S X t X X t X X t X X t X 

— dk = dok = dk = — d). do = "Hq Tik = — Tik "Ho • (112) 
Using the relation [cf. eqs. (97), (111)] 

Xt X X xt 

Hq = Hq Hq = HqHq I4 = I4 (113) 

d 



we find ioT dk A di, {k ^ I, k,l = 1, . . . , 4) 



t X xt 



dk Adi = - dkdi = dki = ^ { di, di - d^ dk 



1 



:t X xt 



— 2 y^k^i " '^i 'Hk ] — — 'Hk'^i — 'Hi'Hk 

= 2 {UkAniY^hi = C2{n)„i'' dM • (ii4) 

Finally, only studying objects with index k = —1 remains to be done. In view 
of eq. (21) it turns out to be useful to define and to study the object 

S X xt X xt X X xt X xt X 

'H[^i) imphcitly defined by means of eq. (115) is a row vector (1x6 matrix). 
For the identical Bogolyubov-Valatin transformation one quickly finds = 
(1, 0, 0, 0, 0, 0). From now on we denote by the symbol T-L the (quadratic) 6x6 ma- 
trix H with row number k = — 1, . . . , 4 and column number m = — 1, . . . , 4 matrix 
elements T-L^k^^ Consequently, for the identical Bogolyubov-Valatin transformation 
we have 

Hi = le. (116) 

For eq. (115), we can also write (emnpgrs is the 6-dimensional totally antisymmetric 
tensor, e(„i)oi234 = 1; we have raised here the indices by means of Ski and not using 
9ki) 

^(-1) = ^{-i) 'Hn'HpHq'H^'Hs ■ (117) 
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Eq. (117) can be transformed [taking into account eq. (99)] to read 

-Tj npqrs 

^"'^ ~ ~5! ^(-1) 

X ^i^'^Hf^^f^Hr^i^'^ ^n'Sp. a,'a,, a,, . (iis) 

7i'.p',g',r',s' — — 1 

Using eq. (22), taking into account eq. (121), and applying the Laplace expansion 
of a determinant [cf. eq. (C.13) in our Appendix C] we find 

= ^ e(_]^) ^'^ T^^" ■''Hg'' ■'H^' ^7/^* ■* e„/p/g/^/^/'" am (119) 



Consistency with eq. (116) now requires that detL = 1, i.e., L E S0{6]R). This 
is closely related to the fact that the unitary group implementing the nonlinear 
Bogolyubov-Valatin transformations [i.e., SU{4)] is simply connected. From the 
eqs. (100), (120) we find in total 

H = L , LL^ = le . (121) 

One can convince oneself that eq. (Ill) also applies for the whole index range 
k,l = -1,...,4. 

In analogy to eq. (110) we can write = zdi(i2^3'^4 = idi A ^2 A ^3 A ^4) 

— c^(-i)o — 2 1 '^(-1) '^0 ~ ^ ^ (-1) 



1 



2 A?/o)'''aM = C2CH)(_i)o'^ aM . (122) 



Furthermore, we have 



^ ^ ^ ^ d ^ d ^ X J 

— id2d^di = —id2 A A ^4 = — 



1 



t X X t 



(^{-1)1 = -j^\d{~i)di — di d{-i) 
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2 AHO'^'eM = C2(H)(_i)/''£m, 



(123) 



idid^di 



^ d ^ d ^ 

idi A ds A d^ 



— C?(-1)C?2 



^^(-1)2 = 2 1 '^(-1) '^2 - ^2 '^(-1) 



- I 'H(„i)'H2 - "Ha 'H(-i) 



2 AHaj'^'aM = C2(H) 



(-1)2 



M 



Cm 



(124) 



-idid2di 



s d s d s 

—idi A d2 A di 



d{-i)d3 



t(-l)3 



2 A?/3)*'aM = C72(?{) 



(-1)3 



M 



Cm 



(125) 



idid2dz 



d d 

idi A d2 A ds 



(i(_i)(i4 



t(-l)4 



^4 — (^4 '^(-1) 



- I 'H(_i)'H4 - 'H(_i) 



Xt X Xt X 



2 {ni-i)An^Y'dM = C2{n) 



(-1)4 



M 



Cm ■ 



(126) 



III. Final result 

Summarizing the above study of two-mode nonlinear Bogolyubov-Valatin transfor- 
mations we can say that these can be studied in terms of nonlinear basis transforma- 
tions of the Clifford algebra C(0, 4). The 15-dimensional space of nontrivial opera- 
tors (the unit operator omitted) in the fermionic Fock space can by spanned in terms 
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of biparavectors. In terms of biparavectors of the Clifford algebra C(0,5), nonlin- 
ear Bogolyubov-Valatin transformations can be described by means of the equation 
[Lm = 'Hm are the row vectors with row number m of the matrix L = 1-1 ^ SO{Q \ M).] 

d M = U Cm = Xm^ Cat 

= 2{nm.,Anm.y £n = 2{L^,ALmy = 02(1) ,f d N ■ (127) 

In a more compact notation we can write 

i) = [jd [J^ = x & , (128) 

where D and C are column vectors with biparavector components (To achieve a 
more compact graphical display, we give here the Hermitian conjugate of C .): 



t 



;t :t xt ;t 



C — ^C(_;^-)Q, C(_;^);^, C(_]^)25 ^(_l)3, C(_]^)4, 

xt ;t ;t ;t xt xt xt xt xt xt 

•^01' '^025 '^03' '^04' '^12' '^13' '^14' '^23' '^24' '^34 

= (^CiC2C3C4, ZC2C3C4, — ZC^CaQ, iCiC2C4, —iCiC2C^, 

Cl, C2, C3, C4, C1C2, C1C3, C1C4, C2C3, C2C4, C3C4j (129) 

and X denotes the 15 x 15-matrix x = C*2(-^)- In component notation, we have 

Xm n = (130) 
Finally, as a special case of eq. (127) we display again eq. (109): 

4 = -Xok^CN = - C2{L)q^^ cn (131) 

expressing the new (transformed) generators of the Clifford algebra C(0,4) in terms 
of the original (untransformed) biparavectors. 

Eq. (127) can be reformulated in a different way. Consider now D and C not 
as column vectors in the 15-dimensional biparavector space /\^{Vq) but rather as 

antisymmetric 6x6 matrices with biparavector entries for which we now write V 

and C , respectively. For example, the antisymmetric matrix V has the matrix 
elements 



1 



t 



^ki = - \ d^,di -di dk] . (132) 
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Eq. (127) then equivalently reads 



V = LCL^ . (133) 

Results corresponding to eq. (133) have been obtained earher by Hristev, [142], ext. 
French version in Rev. Roum. Math. Pures AppL, p. 637, eq. (15.7), [143], p. 176, 
eq. (4.1), ten Kate [102], p. 183, eq. (9), and Buchdahl [94], p. 363, Sec. 8. 

As the matrices L and — L describe the same Bogolyubov-Valatin transforma- 
tion (i.e., the same matrix x) the group of two-mode nonhnear Bogolyubov-Valatin 
transformations is equivalent to the group S'0(6;M)/Z2 (cf. also [102], specifically 
Sec. II and Theorem 1 therein; for related considerations see [138], Vol. 2, Sec. 9.4, 
p. 316, [104], [111], Sec. 5.3, p. 162). 



IV. Further analysis 

On the basis of the equation [cf. the eqs. (105), (131)] 

4 = -Xo/'^M (134) 

and taking into account the eq. (36), we can quickly rederive and reformulate the 
expression for the anticommutator (57). We find 

4 

m=—l 

- « Xok Xoi ^MN Cp (135) 

and, consequently, 

= dkdl + didk= 2 gki = -2 6ki li 

= - 2 Xofc*^ Xoi M I4 - 2i Xo/^Xo/^ ejv/Tv^ cp . (136) 
The eqs. (59)- (63) can be written compactly as 

Xo/^Xo/ M = Ski , (137) 
Xok^Xoi"" ^MNP = 0. (138) 

Eq. (135) now reads 

4 

dJi = - Yl Xor^Xor' im. - Sm I4 (139) 

m=— 1 
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-Xki"" - SmU, (140) 



m=—l 





= 2 Ski 




= 0, 


^MNP 


= 0, 


^1 ^MNP 


= . 



Introducing the complex vectors (e';,)"^ = f^f^^\ i^^k^^\ ^fc°'^\ '^1^'^'°'') 

(1|2) (2|1) (0|1,2) (2|2) (1,2|1) (1|1,2) (1,2|2) (2|1,2) (1,2|1,2)\ _ . . 

Re e'fc = Xo(2fc-i) m, Im = xo(2fc) m) the eqs. (137), (138) can be written as 

(142) 
(143) 
(144) 
(145) 

Consequently, the 15-component complex vectors e'^^, e'g should be two orthogonal 
isotropic vectors of length a/2 fulfilling the additional conditions (144), (145). 

Eq. (137) which is the generalization of the condition (1) characteristic for the 
original (linear) Bogolyubov-Valatin transformations can be understood as an or- 
thonormality condition in the 15-dimensional biparavector space /\^{Vq) and eq. 
(138) is (for k = I) a. decomposability condition. Consequently, any set of four 
operators dk fulfilling the Clifford algebra analogue of the canonical anticommuta- 
tion relations should be a set of four decomposable orthonormal biparavectors. A 
related result has been obtained earlier (in a somewhat more general form, i.e., for 
pentades) by Haantjes, [144], p. 51, stelling 5 [proposition 5], and a corresponding 
comment can also be found in ref. [94], §14, p. 269, above of eq. (14.8). For a further 
discussion of these aspects see Subsec. 6.3. 

The eqs. (137), (138) can be further generalized. Inserting the equation [cf. eq. 
(127)] 

du = Xm^ (146) 
into the anticommutator (38) we find in generalization of the eq. (137) 

Xpip2 Xqiq2 M = ^p\q\^P2q2 ~ ^piq2^P2qi ) (147) 

and also 

X pip2 Xm qiq2 — ^piqi^P2q2 ~ ^piq2^P2qi • (148) 
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The generalizations of eq. (138) read 

Xj'^' Xk"" Xl'' ^MNP = ejKL, (149) 

euNP X*j x\ X^L = ^JKL ■ (150) 

If all the six indices ji, j2, ^i, ^2, ^1, and I2 are chosen pairwise differently, the eqs. 

(149) , (150) can be transformed to read 

detx = 1 • (151) 

This sharpens the condition det x = ±1 that can be derived from the eqs. (147), 
(148). If two of the six indices ji, 72, ^1, ^2, ^1, and I2 are equal, say ki = li = j 
(The cases ri = r2 are trivial and need not to be discussed.), the eqs. (149), (150) 
can be transformed to read [using the orthogonality condition (147), (148)] 

Xjk'^X.i'' eMNP = 0. (152) 

Finally, we study the connection between the conditions (147), (148), and (149), 

(150) . Using the known expression for x in terms of L [eq. (130)], the eq. (147) can 
be written as 

C2(L)C2(L^) = I15. (153) 

Taking into account the compound matrix relation (Laplace expansion of a deter- 
minant; for the notation see Appendix C) 

C2{L)C,{Ly = li5 (154) 
we can immediately conclude that 

C2{L) = C^{L^y . (155) 
This equation can also be written in the following form 

XPQ = l^ ep'''' eQ^^"" XKLXMN (156) 

which can also be derived from the eqs. (149), (150) using the eqs. (147), (148). 



4 Diagonalizing Hamiltonians 

Relying on the insight obtained in the previous section into the structure of nonlinear 
Bogolyubov-Valatin transformations, in this section we will study the diagonaliza- 
tion of certain fermion and spin Hamiltonians. In principle, the methods of this 
section can also be applied, with certain modifications, to other quantities of phys- 
ical interest, for example, the density matrix which is widely being investigated in 
quantum information theory. 
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4.1 Diagonalizing two-fermion Hamiltonians 

Let us now look at an arbitrary two-fermion Hamiltonian H' . 

/,(0|2) 

+ alaioa + /i^^l^'^^ 4«i«2 + /i^^'^l^'^) aja^aiaa (157) 

We assume it to be hermitian [H' = H''^). From the hermiticity condition the 
following 10 relations derive. 



him 


= him 


(158) 


hm 


= him 


(159) 


him 


= him 


(160) 


hi^m 


= - /i(0|l,2) 


(161) 


hiMD 


= him 


(162) 


him 


= him 


(163) 


him) 


= him 


(164) 


himi) 


= - /i(l|l,2) 


(165) 


/,(1,2|2) 


= - /i(2|l,2) 


(166) 


himm 


= himm 


(167) 



Taking into account the equations (158)-(167) (in particular, that h^^^^\ h^'^^'^\ 

/j^(i.2|i,2) real), eq. (157) reads 

- I Re + ^ ci - « Im [h'^^\^^ + ^ /i^^'^l^)^ 62 

-.R.(.-.-i.<-).-Un,(.-.-i..-), 
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- h^'^'^ c,c, - ^ Im (/.(^-^l") + C1C3 

+ 1 Re (/i(1.2|0) _ ;,(l|2)^ g^g^ ^ ^ (;^(l,2|0) ^ ;^(1|2)) g^g^ 

+ ^ Im - /^(il^)) C2C4 - ^ C3C4 

- ^ Re /i^^'^l^) «ciC2C3 - ^ Im /i^^'^l^^ iciC2Ci 

Zj Zt 

- ^ /i(l'2|1.2) iciC2C3C4 (168) 
= |^/,(0|0) + 1 ;,(1|1) + 1 /,{2|2) _ 1 /,(1,2|1,2) j 

-'-tTv^Ird^ . (169) 

Here, the subscript Vq indicates that the trace operation is carried out with respect 
to the six-dimensional paravector space Vg and the antisymmetric matrix Y has the 
exphcit form (To simphfy the display we have omitted the lower triangle matrix 
elements.) 



33 



-Im /i(^'2|2) 

Re(2/i(l|0) + /i^2|2)^ 





Re 

Im (2/i(i|o) + 




Im 

Re [2h^m _ /,(i,2|i)) 
Im(/i(i'2|o) + /iW2)) 

_Re (/i(1.2|0) + /,(1|2)) 





-Re 

Im (2/i(2|o) _ /,{i,2|i)) 
_Re (/,(i-2|o) _ him) 
-Im - /i{i|2)) 

/,(2|2) 



Eq. (169) expresses the Hamiltonian (157) in terms of a biparavector [tryg J; 

plus some constant]. This biparavector stands in an one-to-one correspondence to 
an antisymmetric matrix (exactly as this is the case for ordinary bivectors). The fur- 
ther analysis of the Hamiltonian H' will be based on this correspondence. It should 
be mentioned here that recently Uskov and Rau [145], Appendix B, p. 022331-8 [see 
the first equation below from eq. (B.l)] have given a representation of the Hamilto- 
nian analogous to eq. (168). To explicitly see that both formulations are equivalent 
one must rely on our equations (D.7)-(D.12) given in the Appendix D. 



The matrix Y can be brought to the standard block diagonal form given by the 
matrix Z 



( 

-z^(-i)o 




V 



^-1)0 


















1^12 



-V\2 














\ 







/ 



(171) 



by means of an orthogonal transformation L: Z = LY L . Eq. (169) then reads 



H' 



him 



him 



him 



/,(1,2|1,2) ^ 



- - try, ( ZV 



(172) 



him + i 

2 



1 him _ 1 /,(i,2|i,2) I ^ 

2 4 ' 



+ 2 '^(-l)O " (-1)0 + 2 ^^12 " 12 + - ^^34 " 34 

( him + i + i h^m _ 1 /,(l-2|l,2)^ 1 

V 2 2 4 J ^ 



(173) 



+ 2 '^(-i)o ^c?i44c?4 - - ^12 did2 - - z/34 d^di 

h^'\'^ + \ (/.(^l^) - z/12) + \ {h^'\'^ - ^34) 



(174) 



(/.( 



1,2|1,2) 



2z^(-i)o) 
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+ ^12 b\bi + 1/34 62^2 - 2 z/(„i)o b\blbib2 (175) 

with V = LCL'^ [cf. eq. (133)]. The Hamiltonian is given here in terms of three 
(commuting) Cartan elements: b\hi, ^2^2? ^i^2^i^2 [This is the maximal number for 
the group 5*0(6; M) ^ SU (4).]. For a related consideration see ref. [146]. In a general 
situation, the number of nonvanishing pairs of eigenvalues -^wr [K = (—1)0, 12, 34] 
of the (similar) matrices Y, Z is called the length ([147], vol. 2, Sec. 4.1, Exercise 
26, pp. 52-53), rank ([148], Vol. II, Sec. 9.3.5, p. 331), or mass ([139], p. 67) of the 
related bi(para)vectors. We will use the term rank to diminish the risk of any mis- 
understanding in any physics-related context. We should point out here that due to 
the one-to-one relation between antisymmetric matrices and bivectors on one hand 
and the one-to-one relation between antisymmetric matrices and biparavectors on 
the other hand results available in the literature concerning bivectors carry over to 
biparavectors (we are concerned with) with little change. The standard orthogonal 
decomposition of bivectors^^ relies on the fact that any antisymmetric matrix can 
be brought to block diagonal form [the diagonal blocks are proportional to J ); 
see, for example, [151], Chap. XI, §4, various pp. in the different Russian editions 
(Vol. 2, pp. 12-18 of the English transl.)]. For the general case of three mutually 
different pairs of eigenvalues ±zz/(_i)o, ±ii^i2, ±^1^34, of the (similar) matrices Y, 
Z a geometric algebra formalism to calculate these eigenvalues and the orthogonal 
transformation L to transform the matrix Y into the matrix Z exists (We will not 
review it here. Cf. [152], Chap. 3, Sec. 4, p. 78, [153].). This general situation, for 
example, will apply for Hamiltonians of (quasi)fermions which have been obtained 
from Jordan- Wigner transformations of two-spin- ^ systems (Those will be discussed 
in some detail further below.). However, in the present subsection primarily we have 
in mind the interpretation of H' as the Hamiltonian of two (equivalent) fermionic 
modes. Consequently, we expect for a physical Hamiltonian that two of the eigen- 
value pairs [with an eye to eq. (175) chosen to be ±^12, ±^^'34] should be equal 
zi/12 = W34 [i.e., the Hamiltonian should not change under exchange of the fermion 
indices 1 and 2 - this entails further identities for the coefficients in eq. (157)]. Con- 
sequently, eq. (175) exhibits a residual 0(2) ® 0(2) symmetry related to the twofold 
degenerate eigenspaces of the matrices Y, Z related to the eigenvalues 11^12 = ivm 
and —ivi2 = —ii^sA, and a certain arbitrariness exists in constructing the orthogonal 
transformation L (and the related canonical transformation) [56], Sec. 17.3, p. 222 
(In this reference, this fact is expressed in terms of a bivector expansion.). 

The eigenvalues A = ^iux [K = (—1)0, 12, 34] of the matrix Y can be calculated 
by means of the characteristic equation (For the trace operation we omit here the 
subscript Vg because no misunderstanding may occur.) 

= det{Y-\le) 

i6[136], all eds., Chap. II, §11, pp. 35-36, [149], [148], Vol. II, Sec. 9.3.5, p. 331, [139], [150], Sec. 
1, p. 184, [134], p. 103 (p. 2251 of the English transl.). 
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= - -(trF^) + -{(tiY^Y - 2trHA2 + detF (176) 
2 8 L J 

which is a cubic equation in terms of A^. The invariants of the matrices Y, Z entering 
the above equation are 

tiY^ = tiZ^ = -2 (z/(%)o + ^^12 + 4) , (177) 

tr = tr = 2 (z/f.^^o + u^, + ul,) , (178) 

detF = detZ = (i^{_i)oz^i2i^34)^ • (179) 

By calculating the discriminant A of the equation (176) (We do not display here the 
explicit expression of it, see any standard reference on cubic equations.) one can 
check if indeed two of the eigenvalue pairs A = ±iiyK agree as expected. If this is 
the case the eqs. (177)-(179) read 

tr Y^ = tiZ^ = -2 + 2z/2) , (180) 

try^ = trZ^ = 2 + 2i/^) , (181) 

detF = detZ = (i^(_i)oi^^)^ , (182) 

and one can quickly find from these equations the eigenvalues A = ±zz/(_i)o, A = 
±1/ = ±ii/i2 = ±i^'34 without the need to resort to the standard machinery of solving 
a cubic equation in general. From the eqs. (180), (181) one finds for the squared 
eigenvalues A^ 



n-i)o = ^ -t^^' T 2^3trF^-i(trF2)2l . (134) 



Note, that 3 tr F"^ — ^ (tr F^)^ = 4 (^u"^ — ^f_i)Q^ ■ Consequently, in the eqs. 
(183), (184) on the r.h.s. the upper sign applies to weak (quasi)fermion coupling 
(z/2 > z^(2_i)q) and the lower sign to strong (quasi)fermion coupling (i/^ < iy'^_^-jQ). 
Which sign applies can be determined by calculating the determinant (or the Pfaf- 
fian) of the matrix Y. 



One may wonder why the analysis of the 4- level Hamiltonian H' [eq. (157)] which 
can be represented by means of a 4 x 4 matrix leads to a cubic equation and not to 
a quartic one. The answer is as follows: Writing the Hamiltonian H' in terms of a 
biparavector plus some constant [eq. (168)] amounts to splitting the Hamiltonian H' 
into a traceless part plus some diagonal term (proportional to I4). Consequently, as 
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this way one degree of freedom for the eigenvalues of the Hamiltonian H' has been 
separated out, the remaining characteristic equation for the traceless part of H' is 
reduced in order by one degree to a cubic equation. 

Finally, let us dwell on certain general considerations. Qualitatively, the class 
of physical Hamiltonians H' (with 1^12 = z/34 = v) can have rank 2 (z^{-i)o = 0) or 
3 (z^(-i)o 7^ 0) - interpreted in terms of the biparavector(s) related to it. For the 
generic case ^^(-1)0 7^ (i.e., rank 3) it is clear that it is impossible to write the 
Hamiltonian H' as the sum of two noninteracting quasiparticle fermion oscillators, 
i.e., it is impossible to define canonical transformations of the fermion creation (aj, 
Sg) and annihilation (Si, 0,2) operators such a way that the Hamiltonian H' can be 
written in terms of noninteracting quasiparticles [(quasi)fermions] (61, 62)- However, 
if 

detF = (PfY)^ = detZ = (Pf Z)^ = (z/(_i)o z/^)^ = (185) 

(i.e., z^(-i)o = 0) this is possible. 

We conclude this subsection with some comments on the related literature. The 
Hamiltonian (175) with z/12 = 1^34 = v has been studied for fermion systems in 
[48, 50, 52, 154-158], while in the references [159-161] the somewhat more general 
situation where not necessarily 1^12 = ^'34 is being considered. Particular versions of 
the general Hamiltonian (157) have been studied in [154, 155, 162]. However, the 
case studied in [162] can be treated by means of linear Bogolyubov-Valatin trans- 
formations. 

4.2 Nonlinear spin transformations 

— Diagonalizing spin Hamiltonians 

4.2.1 Two-spin-| systems 

It is well-known that fermion systems can be related to spin systems by means 
of Jordan- Wigner transformations [163]. This, of course, also applies to the two- 
fermion system under consideration. From a group-theoretical point of view, each 
elementary spin (|) corresponds to some SU{2) subgroup of some larger group. 
Inasmuch as spin operators (Pauli operators) of different spins commute among each 
other a system of n spins corresponds to n pairwise commuting SU{2) subgroups 
of the larger group related to the system under consideration. In our setting, the 
transition from fermion operators to spin (Pauli) operators is related to the choice of 
two commuting su(2) Lie subalgebras within the su(4) Lie algebra of the two- fermion 
system. By means of this choice a two-fermion system is mapped to a two-spin- 1 
system. Such systems have been studied from various points of view in the past 
[164-172], [2], Chap. 2, Sec. 13, pp. 52-56, Table 8, pp. 294-299, [173], [174], Chap. 
5, p. 120, [84, 175-186]. We leave here aside the broad range of studies performed 
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in recent years concerning the problem of entanglement in two-spin- 1 systems - 
the "harmonic oscillator of quantum information theory". It is clear that there are 
many possible Jordan- Wigner transformations (This question seems not to have been 
studied systematically in the literature so far.). Here, we will follow [84] with some 
modifications arising from certain esthetic considerations related to the equations 
(24)-(27). We want to choose a version of the Jordan- Wigner transformation that 
is fairly symmetric with respect to the mode number indices involved. The authors 
of ref. [84] define (p. 13835, above of eq. (5); [S^, Si] = iSl = SI ± iSf) 



St 



I4 — (1 + i) 0,20-2 



(186) 



Si 

SI 

St 



Oi 



'1 - i) ala2 



oloi - 2 ^4 



I4 — (1 — i) ojoi 



(187) 
(188) 
(189) 



S' 



02 



I4 — (1 + i) ojoi 



0302 



1 

- I4 . 

2 



(190) 



(191) 



One can convince oneself without any difficulty that the inverse Jordan- Wigner 
transformation from the spin-^ (Pauli) operators to fermion operators is given by 
the following equations. 



Oi 



02 



i S^ 
I St 
i S2 
i St 



^(1 + I4 



2(1-^) I4 



;i - z) SI 



:i + s^ 



[l-l) l4 + (l+^) ^1 

(1 + z) I4 + (1 - t) S'^ 



(192) 

(193) 
(194) 
(195) 



In difference to ref. [84], we define spin-| operators by means of the following equa- 
tions [Our choice is related to the choice for the eqs. (24)- (27).]. 



c+ _ " 



[l + i) I4 — 2 0,302 



[1 — i) I4 — 2 0302 



(196) 
(197) 
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St 



2 



V2 

i 



a\ai 



0-2 



1. 



1 — z) I4 — 2 a\ai 
{l + i) I4 - 2 a[ai 



And the inverse transformation is given by 

1 



a2 



V2 
1 

71 
1 

71 
1 

72 



5+ [I4 - 2z , 

5f [I4 + 2z 51] , 

^2+ [I4 + 2i SI] , 

[I4 - 2z . 



(198) 
(199) 

(200) 

(201) 

(202) 
(203) 
(204) 
(205) 



It is clear that any two-spin- 1 Hamiltonian can be transformed by means of the 
eqs. (186)-(191) to the general form (157) of the two-fermion Hamiltonian. It is 
furthermore possible, as discussed above, to bring any two-fermion Hamiltonian 
to the special form (175). Of course, also the quasifermion operators 61, 62 can be 
related to quasi-spin-| operators T^, T^, {k = 1, 2) by means of a Jordan- Wigner 
transformation. In analogy to the above equations we can write 



bi 



--j= T+ [U4 - 2z T^] , 



1 

72 
1 

71 
1 

7! 



Tf [I4 + 22 T^] , 
T2+ [I4 + 2z Tf ] , 
[li-2iT^] . 



(206) 
(207) 
(208) 
(209) 



Then, the quasifermion Hamiltonian H' [cf. eq. (175)] can be written as a two-quasi- 
spin-i Hamiltonian. It reads 



H' 



/,(0|0) + 1 _^ 1 f^(2\2) _ 1 ;^(1,2|1,2) | 
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- (i^(-i)o + 1^12) n - (^(-1)0 + ^^34) T| + 2 z/(_i)o Ti^Ta" . (210) 

Consequently, any two-spin- 1 Hamiltonian can be brought to the above form. A 
related result has been found (for a somewhat restricted class of two-spin- ^ Hamil- 
tonians) in [2], Chap. 2, Sec. 13, pp. 52-56. The unitary transformations reducing 
any two-spin- 1 Hamiltonian to the form (210) correspond, in general, to nonlinear 
spin transformations. Such transformations have been considered in a somewhat 
different context in [187], see p. 1186, eq. (6a), [188-193]. 



4.2.2 Single-spin-| systems 

Two-fermion Hamiltonians are not only related to two-spin- 1 systems but can also 
be related to a single-spin- 1 system (A concrete choice for the latter relation can be 
found in [194].). Consequently, any system which is given in terms of one of these 
representations can also equivalently be formulated in terms of the other two. The 
study of single spin-| systems has a long history (For certain theoretical aspects 
see [195-198] and references therein.) and they have found recent attention in the 
field of quantum computation [199-202]. To make the relation of any two-fermion 
system to a single-spin- 1 system explicit we will follow here [194]. We will apply the 
same line of reasoning as for a two-spin- ^ system applied above. The three (radial) 
spin-| spin operators J"*", /~, P can be given in terms of the fermion creation and 
annihilation operators the following way ([194], p. L506, eqs. (19), (20)). 

/+ = v/3 02 + 2 a\ai (211) 

/- = 4 + 2 a{a2 (212) 

P = = 2alai + 4«2 - ^ I4 (213) 



The inverse transformation reads ([194], p. L506, eq. (21)): 

1 

Vs V2 

1 



^ , (214) 



di = — ^rrr , (215) 



^ + , (216) 



a. ^ ^-U + F) /-. (217) 

Again, any single-spin-| Hamiltonian can be transformed by means of the eqs. (211)- 
(213) to the general form (157) of the two-fermion Hamiltonian. It is furthermore 
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possible, as discussed above, to bring any two-fermion Hamiltonian to the special 
form (175). Then, also the quasifermion operators bi, 62 can be related to quasi-spin- 
I operators J+, J~, by means of a Jordan- Wigner transformation. In analogy 
to the above equations we can write 

J+ = 73 62 + 2 blk , (218) 

J- = V3bl + 2 b\b2 , (219) 

.r = l[J^,J^] = 2blk + blh - ^14, (220) 



and 



b\ = - ^ J+.rJ+ , (221) 
61 = - ^ J~J'J- , (222) 



v^^' V2' 



b\ = ^ J+ ( + r] , (223) 



1 /I 



62 = ^ (^-14 + J^J (224) 

The quasifermion Hamiltonian H' [cf. eq. (175)] can be written as a single-quasi- 
spin-| Hamiltonian. It then reads 



2 2 4 

+ ^ (17 z/(_i)o - 5 1/12 + 22 1/34) + ^ '^(-i)o (^')' 

+ ^ (/^(-i)o - + 2 1/34) (J')' . (225) 

Finally, we would like to mention that on the basis of the two above (generalized) 
Jordan- Wigner transformations for two-spin- 1 and single-spin- 1 systems relations 
can be established between spin-i and spin-| operators^''. Combined with the full 
range of nonlinear Bogolyubov-Valatin transformations discussed we thus obtain a 
large manifold of expressing single-spin- 1 operators in terms of two-spin-| operators 
and vice versa. 



^^There seems to have been done some related research in the past [203]. 
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5 SU{A), 50(6; M), 50(6; M)/Z2 transformations and 
their parametric relations 

In the previous sections unitary SU{4) transformations U, orthogonal 5*0(6; M) 
transformations L, and 5*0(6; M)/Z2 transformations x have played an important 
role. So far, we have not discussed (except for the expression of x in terms of L) their 
concrete mutual relationship which is of importance for any explicit calculation. In 
this section we finally will consider this technical problem. To set the frame for this 
discussion we will first specify which sort of parametrization for the matrices U and 
L we are going to use. 

Let us write for the unitary 4x4 matrix U = U ({A}; {c}) implementing a given 
nonlinear Bogolyubov-Valatin transformation (Tq is a complex number while T'"^™^ 
are the matrix elements of a complex antisymmetric 6x6 matrix T.) 

U = To I4 + T^' Im = ToU + ^ T^'"-"' lm,m, . (226) 

From the unitarity condition UU^ = I4 follow the equations 

|To|2 + T^^ Tm = 1 , (227) 

— Tq Tp + Tq Tp — Tp^"* Tmp2 + T^p^ + ^ T^'^ T eMNP = . (228) 

These equations have been given earlier (in some less general notation) in [204], p. 
243, eq. (35). For another version of these equations see eq. (332). For any given 
unitary matrix U the coefficients Tq, Tm can be calculated by means of the equations 
[cf. eq. (36)] 

To = ^ tr U , (229) 
Tm = -7 tr (Sm u) . (230) 



4 

For the orthogonal 6x6 matrix L we rely on the Cayley representation 

Iq + A 2A 2 , , 

L = — 7 = 16 + 7 = -le + 7 231 

in terms of a real antisymmetric matrix A (We disregard here all problems of the 
Cayley representation related to any eigenvalues = -1 of the orthogonal matrix L.). 
The matrix L can be expressed as a sum over a finite number of powers of the matrix 
A. For an explicit expression see eq. (274) further down. In turn, the antisymmetric 
matrix A can be expressed in terms of L as 

L-te _ T(-) 
" LTTe ~ l6 + LW ' ^232) 
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= ^{L + L^) , (233) 



L(-) = ^{L-L^) . (234) 

The antisymmetric matrix A can also be expressed as a sum over a finite number 
of powers of the matrix L, however, already for the present case these expressions 
are quite involved and, therefore, we will not display them here (for details see [205]). 

To simplify navigation through the present section we now give in Figure 1 a 
schematic overview of its content. Each arrow in Figure 1 stands for an equation 
expressing the quantity at the end point of the arrow by another one at its starting 
point. 




Figure 1: Schematic overview over the content of Sec. 5. 
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List of reference points for Figure 1: 

Arrow 1: Eqs. (229), (230). 

Arrow 2: Eq. (226). 

Arrow 3: Eq. (232) and ref. [205]. 

Arrow 4: Eq. (231), Subsec. 5.2.1, eq. (274). 

Arrow 5: Subsec. 5.1, eqs. (248), (249). 

Arrow 6: Subsec. 5.2.2, eq. (281). 

Arrow 7: Subsec. 5.2.1, eq. (268). 

Arrow 8: Subsec. 5.1, eq. (250). 

Arrow 9: No equation displayed, can be composed of the equations corresponding 
to tlie arrows 6 and 1 (or tlie arrows 3 and 12). 

Arrow 10: No equation displayed, can be composed of the equations corresponding 
to the arrows 5 and 3 (or the arrows 1 and 11). 



Arrow 


11: 


Subsec 


Arrow 


12: 


Subsec 


Arrow 


13: 


Subsec 


Arrow 


14: 


Subsec 



5.1 From U G SU{4) to L e 50(6; R) 

In this and the following subsection, we will study the explicit relation between the 
5*0(6; M) transformations L and the corresponding unitary SU{4) transformations 
U ({A}; {c}). In the present subsection, we will assume that U = U ({A}; {c}) is 
known and derive from it an expression for L. We will perform this task by relying 
on the discussion presented in [206, 207] (for follow-up work see [208]). 

To begin with, let us start by repeating some elements of the very lucid account 
given in [207] (We will follow here also the notation used in this article with slight 

modifications.). Consider a set of six linearly independent antisymmetric (complex) 

+ + J. 

4x4 matrices^^: = — , k = —1, . . . , 4. These matrices represent spin (space) 
bivectors (e A^(*^4)' where C4 is the spin space) and can be chosen to obey (e"^'='^, 



^^For an early discussion of the role of antisymmetric matrices in the present context see [209], 
related discussions can be found in [210, 211], [212], Chap. IV, Sees. 5-9, pp. 36-40, [213] (There are 
many other papers by Schouten and collaborators containing related but less focused material.), 
[144, 214] (very clear mathematical accounts), [215], § 11, p. 49-52, [217], [94] (also note [103]), 
[104], [218], JIeK^. [Lekts.] 13, pp. 258-299 (pages 247-285 of the English transL), and in [219-222]. 
A somewhat related discussion for the case SL{2;C), 50(1, 3; M) is given in [223]. A different 
approach of dealing with the homomorphism between the groups 5C/(4) and 5*0(6; R) related to 
the physical problem of the 3-particle problem is discussed in [224], Sees. VI. A, VLB, pp. 557-559. 
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eabcd, e^^^"^ = ei234 = 1 is the completely antisymmetric tensor operating in spin 
space C4) 

([207], p. 10, eq. (1), an example for these matrices can be found in [206], p. 814, eq. 
(12); also see our Appendix D). One can then define antisymmetric 4x4 matrices 

" ~ T + - 

= — (the spin space Hodge duals of by writing 

ab 1 

abed 



r.) = -2^"^"' [^k)^^ ■ (236) 



Then, eq. (235) can be written compactly as 

^tr(f,r,) = 6ki. (237) 

+ 

From this equation one recognizes^° that the matrices , F^ obey the equation^^ 

f^r, + f,F, = F,f, + F,f, = 24/14. (238) 
One can also derive the following useful relations^^ 

^k)^S^'T = '^^'^'^ ~ ^^^^^ 

- \ ab ^- ^ cd 

) f fM = 2 e"'""^ , (240) 
(f'') = 2eabcd. (241) 

ab \ / cd 



^^Incidentally, the minus sign on top of the F symbol should not be confused with the sign for 
complex conjugation ~ a longer bar. 

20Cf. [210], p. 410, eq. (2.14), [212], p. 36 (p. 216 of the whole volume), bottom of the page, 

[213], p. 184, eq. (60), [144], p. 48, eq. (2), [214], p. 139, eq. (1.5), 

[94], p. 368, eq. (13.8), [207], pp. 10-11, [221], Subsec. 2.1.1, p. 19, eq. (49), 

[222], p. 11 (p. 9 of the English transL), eq. (2), [225], p. 24, Appendix A, eqs. (A2). 

^^Incidentally, such a relation has also emerged in [226], p. 1066, eq. (6). For a somewhat related 
discussion also see [64], in particular, p. 3, eq. (5). 

22Cf. [211], p. 510, eq. (4.10), [213], p. 177, eq. (7), [144], p. 48, eq. (6), [214], p. 138, eqs. (1.2), 
[215], p. 50, eq. (33) (incidentally, in the preceding eq. (32) 774^ and 77^5 should correctly read rjf^ 
and 77^3, respectively), [216], p. 159, eq. (21), [94], p. 367, eq. (13.6), 

[62], p. 93, eq. (2.10), [207], p. 11, eqs. (11), (12), [220], p. 44, eq. (2), 

[221], Subsec. 2.1.1, pp. 17-20, [222], p. 11 (p. 9 of the English transL), eq. (1), 
[225], p. 24, Appendix A, eqs. (A4)-(A6). 
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It is now possible to represent the paravectors in terms of the antisymmetric 



matrices , . One finds^^ 



= -FoF, . (242) 



The choice of Fq singles out a certain spin bivector and in the spin bivector space 
/\^(C4) a certain direction (vector, line). Haantjes calls this direction (line) the 

axis of the pentade c^, k = —1,1..., 4 ([144], p. 51, stelling 4 [proposition 4]). 

+ 

Furthermore, it is possible to choose the matrices F^ , F^ such a way that 

ffc^ = r, (243) 

applies. The condition (243) determines in the (complex six-dimensional) spin bivec- 
tor space A^('^4) ^ (real six-dimensional) subspace which we denote in view of eq. 
(237) by Ar5('^4) — [aAso see the comment below from eq. (247)]. As a techni- 

cal comment, we would like to mention here that F^ on the l.h.s. of eq. (243) is 

an object with upper (spin) indices related by hermitian conjugation to the object 
+ 

F^„ with lower (spin) indices. To cut a long story short, this is related to the fact 
that for bivectors in Ar6(^4) — -'^e? we can raise and lower indices in spin space by 
means of Sab' = ^aWi s""^' = supplemented by an operation of hermitian conju- 
gation [cf. [221], Subsec. 2.2.3, p. 40, table 1, and Subsec. 2.2.2, p. 31/32, TeopeMa 
[Teorema] / [Theorem] 2, eq. (121); [222], p. 15 (p. 12 of the English transl.), table, 
and p. 13 (p. 11 of the English transl.), Teopeivia [Teorema] /[Theorem] 2, eq. (10)]^^. 

By virtue of eq. (236), eq. (243) leads to the condition ( — denotes complex 
conjugation) 

m = (r,)°' = -1.-- (f,)^ (244) 

+ 

which entails that the real part of the matrix F^ should be a selfdual matrix while 

the imaginary part should be an anti-selfdual matrix. We give in Appendix D an 

+ 

explicit example for the set of matrices F^ that obey eq. (243) [Up to the choice of 
signs this set of matrices agrees with that given in [206], p. 814, eq. (12).]. Taking 
into account eq. (243) we find 

t = -f,Fo . (245) 



23 Also cf. [210], p. 410, cq. (2.13), [212], p. 37 (p. 217 of the whole volume), [144], p. 49, eq. 

(7). The analogous equation for an arbitrary biparavector c m [see eq. (336), Subsec. 6.3] is given in 
[94], p. 367, eq. (13.7), [104], p. 2244, eq. (9), [207], p. 11, eq. (18), 

[220], p. 44, eq. (5a), [111], Sec. 5.3, p. 167, eq. (5.21), [221], Subsec. 6.1.1, p. 97, eq. (472), 
[222], p. 12 (p. 10 of the English transl.), the first eq. on the page. 

^''^We are indebted to K. V. Andreev for pointing out to us these subtleties and refer the interested 
reader to his publications [220-222] for further details. 
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S it 

Equations analogous to the eqs. (242), (245) also hold for d k, d with the matrices 

Tf^ replaced by some other matrices F^. These two sets of matrices are related by 
the equation^^ 

f,' = L,iT, = uf,U^, (246) 

= = u r, ut . (247) 

In general, here the matrix L may be an element of the group 50(6; C), however, by 
virtue of the condition (243) the matrix L must be real in our case and, therefore, 
belongs to the group 5*0(6; M). It seems worth mentioning that eq. (244) is invariant 
under the transformations (246), (247) (To see this one has to rely on the Laplace 
expansion of the determinant of the 5?7(4) matrix U.). The eqs. (246), (247) pro- 
vide us with additional motivation to denote the subspace of the spin bivector space 

/\^(C4) we are using by the symbol Me because its invariance group is 50(6; M). 

+ 

The antisymmetric matrices can be understood as the orthonormal basis in the 
spin bivector space /\^(C4) [cf. eq. (237)]. 

From eq. (247) one finds the relation^^ 

Lki = ^tr(F;f,) = itr(F,Utf,u) (248) 

= ^ tr (r,: ) = ^ tr (f , F, U) . (249) 

Inserting now eq. (226) into it and calculating the occurring traces by means of the 
eqs. (E.1)-(E.6) we obtain [cf. [206], p. 822, eq. (53)] 

L = (tI + T^^Tm) le - 2T0T + 2T' - 2iPc\T) 

= {T^ + T^' Tm) Ig - 2 To T + 2 + 2z Pc'^{T) . (250) 

Here, we have introduced the matrix Pc^{T) for the antisymmetric matrix T by 
defining 

Pc\T),i = J^PfT = -i T^^T^eMiVfcz . (251) 

We call the matrix Pc^{T) the (second) supplementary Pfaffian compound matrix. 
We have defined it in analogy to the definition of supplementary compound matrices 

25[206], p. 815, eq. (10), p. 821, eq. (47), [207], p. 12, below from cq. (23), p. 13, eq. (30). Note 
that Stepanovskii has studied the more general case S'L(4;C)/Z2 ~ 50(6; C) and that we have 
interchanged U and compared with the use in [206, 207]. 

26Cf. [206], p. 822, eq. (52), also see [227], p. 766, the equation below from eq. (2.4), [60], p. 13, 
eq. (4.9), [104], p. 2244, eq. (6), [221], Subsec. 2.2.1, p. 22, eq. (62), and [222], p. 12 (p. 10 of the 
English transl.), eq. (4). A related result can be found in [100], p. 139, eq. (3.13). 
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from a determinant [cf. Appendix C, eq. (C.9)]. The matrix Pc'^{T) obeys the 
following equations. 



Pc^{T)T = T Pc'^iT) = PfTUe 



(252) 



[Pc\T)\ 



(tr T'Y-ti 



rp2 y2 _ y4 (^253) 



Under the assumption Pf T 7^ we can derive from the above two equations the 
relation [by multiplying eq. (253) by T] 



Pc\T) 



T 



Pf T 4 



^(tr T^Y-tr T' 



1q - -T"^ tiT"^ + . (254) 



Eq. (252) can be viewed as a special form of the Pfaffian analogue of the Laplace 
expansion of a determinant^^. 



It is clear that the imaginary part of the expressions on the r.h.s. of eq. (250) 
must vanish, consequently the following relations apply. 

Im (T^ + T^Ta/) le + 2 Im (T^) = (255) 

Im(ToT) - Re(Pc2(r)) = (256) 

Eq. (255) originates from the symmetric part of eq. (250) while eq. (256) is derived 
from its antisymmetric part. From eq. (255) we can obtain the equation 

lm{T^) = Im(trr2) . (257) 



5.2 From 50(6; R) back to its double cover S'C/(4) 

Having obtained eq. (250), we will now reverse reasoning and study the relation be- 
tween the 5*0(6; M) transformations L and the corresponding unitary SU{4) trans- 
formations U ({A}; {c}) in the opposite direction. We will now assume that L is 
known and derive from it an expression for U = U ({A}; {c}). In the literature, one 
finds two approaches to this task. One is due to Fedorov and collaborators [235, 236] 
and another one has been given by Klotz [104] (relying on earlier work by Macfar- 
lane [100, 237]). In the first part of the following discussion (Subsec. 5.2.1) we will 
rely to a large extent on the former while in a second part (Subsec. 5.2.2) we will 
describe the approach by Klotz. 

^^Cf. [228], Part 1, Chap. 1, Sec. 3, p. 6, eq. (1.5). This expansion is originahy due to Tanner 
[229], has also been noted later by Baker [230], and proved in [231-233]. For a recent discussion 
see [234]. 
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5.2.1 Plane orthogonal transformations — The approach by Fedorov and 
collaborators 



To begin with, it turns out to be useful to recall that any orthogonal transformation 
5*0(6; M) can be represented in terms of 3 commuting orthogonal transformations 
within mutually orthogonal planes [238, 239], [135], Act. Sci. Ind. 643, §45, pp. 
46/47 (p. 36 of the English transl.). In the Cayley representation of an orthogonal 
transformation L [eq. (231)] each of these plane orthogonal transformations can be 
given in terms of a decomposable (real antisymmetric) matrix A. For such a plane 
orthogonal transformation [obeying = ^ (tr A^) A] the orthogonal matrix L can 
be written as (a 7^ 1; [235], p. 1034, p. 1350 of the English transl.) 

2A 1 

L = le + (te + A) , a = - tr A^ . (258) 

1 — a 2 

Now, comparing the antisymmetric parts of the equations (250) and (258) we im- 
mediately find the relation [The last term on the r.h.s. of eq. (250) vanishes for a 
plane orthogonal transformation.] 

= - To T (259) 



where 

To = ij^l + ^trT^y = ±(l-a)-5 = ± (^1 - 1 tr ' (260) 

in view of eq. (227) [To and T are real for a plane orthogonal transformation.]. For 
a plane orthogonal transformation parametrized by the antisymmetric matrix A it 
holds 

1 - ^ tr = det (le - A) . (261) 
Consequently, eq. (259) can be written as 

T = T(l--trAA ' A = T , ^ , (262) 

V 2 ; + Vdet (l6 - A) ^ ' 

A = - ^ = - sign To (^1 + ^ tr T^^ ' T . (263) 

Finally, any unitary transformation U standing in correspondence to a plane orthog- 
onal transformation L = L{A) can be written as 

U = U(A) = ±[det(U6-^)]"^ (l4-^^£7v) • (264) 
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Having obtained this representation we can now go over to the general case. 
Any real antisymmetric matrix A representing a general orthogonal transforma- 
tion 5*0(6; M) can be written as the sum of three decomposable real antisymmetric 
matrices Ak that represent three (commuting) plane orthogonal transformations in 
mutually orthogonal planes: 



A 



Ai + A2 



A, 



We have 



L(A) 



, Ak Ai = ioi k^l . 
L(Ai) L{A2) L{A^) 



and, consequently, 



U = \}{A) = U(Ai) U(A2) U(A3 



(265) 
(266) 
(267) 



Inserting eq. (264) into eq. (267) and performing the algebra [by means of eq. (36) 
and taking into account the decomposability of the matrices A^ (no summation over 
k here): A^^A^eMNP = ] we find 



U 



\J{A) = To I4 + &p 
= ± [det (le - A)]-^ 

X \{l-iPiA)li- [A^ -i Pc\A)^)lp 
To = ± [det (le - A)]"^ (1 _ ^ PM) , 
Tp = T[det {le - A)]~^Ap-zPc\A)p) . 



(268) 

(269) 
(270) 



Eq. (268) agrees with eq. (27) of [236] (p. 988, p. 258 of the Enghsh transl.). For a 
related result see [104], p. 2245, eq. (13) [or its discussion in our Subsec. 5.2.2, eq. 
(281)]. It also yields an alternative to eq. (250) by relying on the Cayley represen- 
tation [eq. (231)] of the matrix L [cf. eq. (26) of [236] (p. 988, p. 258 of the English 
transl.)] 



A 



ReT 
ReTn 



(271) 



The eqs. (269), (270) can now be inserted into eq. (250) to obtain the following 
further version of eq. (231). 



If 



2A{le + A) 
A) 



det ( If 



(l-^tr A''^ le + A'-[Pc'iA)]'- 



(272) 



Besides more elementary ones we have made use of the following result in deriving 
eq. (272). 



eKMN Re T*^ Im 



det (U 



6 



A) 



1 



Ak tr A^ 



K 



(273) 
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Eq. (272) can be further transformed to read 
L = 1q + 2A 

2(l6 + A) 



det ( If 



A) 



(det A) 1q - A' 



1 - - tr 1 Ifi + 



This form agrees with the general eq. (3) in ref. [235]. 



. (274) 



5.2.2 The approach by Klotz 

The approach by Klotz [104]^*^ (adapted to our case) for relating 5*0(6; M) transfor- 
mations L and the corresponding unitary SU{4:) transformations U ({A}; {c}) starts 
by considering the relation 



cn 



■N 



6X' - 4 5:^5, 



(275) 



which goes back to Pauli . Here, (cat I denotes the matrix elements of the oper- 



ator (4 X 4) matrix c n. Eq. (275) can easily be derived by means of the eqs. (23), 
(242), (245), (239)-(241). Multiplying eq. (275) by U^,'^ andU';, [i.e., certain matrix 
elements of the unitary operator U ({A}; {c}) G S'f/(4)] and performing the sum over 
the indices a, b, c, one can derive the following equation (in operator notation). 



U tr U 



U c/v c 



N 



(276) 



Recognizing, that U c a? can be written as Xn *^cm and taking into account the 
relation tr U = 4 Tq [cf. eq. (229)] we find the expression 



U 



Tn 



16 |To| 



I4 - XNM^^'d'' 



(277) 



By taking the trace of eq. (277) one can derive the relation (tr x = Xm 



1 

16 



[ 1 + tr X 



1 

16 



1 (tr lY - - tr 

2 ^ ^ 2 



(278) 



^^Note, that an analogous approach has been used earher by Macfarlane (leaving aside here the 
problem of the signature in the spaces used) in the cases of 50(4; R) [237] and 5*0(5; R) [100], p. 
145. 

29[240], p. 32 (p. 725 of Vol. 2 of the reprint [241]), eq. (I), 
[242], p. 118 (p. 762 of Vol. 2 of the reprint [241]), eq. (II) Also cf. e.g., 
[243], p. 210, eq. (233), [92], p. 35, eq. (11), [100], p. 137, eq. (2.22), 

[94], p. 361, eq. (4.2), [104], p. 2245, Sec. IV, [96], p. 208, eq. (50), 

[111], Sec. 5.3, p. 167, eq. (5.22), [220], p. 44, eq. (5a), [221], Subsec. 2.1.1, p. 18, eq. (47), 
[222], p. 12 (p. 10 of the English transL), eq. (2'). 
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which, however, still leaves the (complex) phase of Tq undetermined. One can con- 
vince oneself that this result agrees with the corresponding expression derived from 
eq. (269). To see this explicitly one has to rely on the relation 



(tr L)^ - (tr Lf tr 



tr L tr L-^ 



+ - (trL=^)^ - trL^ = 2 (tr L)^ - 2 tr (279) 

which can be derived from eq. (155) by means of taking appropriate trace operations 
on both sides. Further progress in discussing eq. (277) can now be made by relying 
on eq. (36). Finally, eq. (277) is found to read [tr L = tr L^+\ Pc^{L) = Pc^{L^-^)] 



U 



= To 
This entails 



1a - 



16 |To 
1 

8^ 



- [(tr L) le - 2 L^] - 2z Pc\L)}^^d 

2 Tn 



M 



(280) 
(281) 



-M 



1 + i (tr ly - i tr L2 



X 



{L^~^ [(tr L)1g-2 lW] - 2i Pc\L)} 
which stands in correspondence to eq. (270)'^*^. 



M 



(282) 



It now remains to calculate the (complex) phase of Tq. Klotz [104] incorrectly 
states on p. 2245 [above of eq. (13)] that the (complex) phase of Tq can assume the 
(fixed!) values 0, |,7r, ^ only. This statement does not follow (as Klotz suggests) 
from the condition det U = 1 and is also in contradiction to eq. (269). Following the 
reasoning of Klotz [104], Laufer [111] has also discussed (Sec. 5.3, p. 162-179) the 
homomorphism between SU{2,2) and 5*0(2, 4; M) (The difference in signature to 
our case 5'f/(4), 5*0(6; M), can be disregarded for the present purpose.) and makes 
the same erroneous statement"'^. As a matter of fact, Tq can be calculated from the 
condition det U = 1 as some fourth (complex!) root and any root will differ from 
the other three by one of the phases |, vr, ^ (mod 27r) but can be complex itself, in 

'^'^Infinitesimally L^~^ ~ 2A, and for this domain one immediately recognizes that eq. (282) 
agrees with eq. (270). However, we have not attempted to exphcitly check the equivalence of the 
eqs. (270) and (282) in general. 

''^Bottom of p. 169; incidentally, the problem of the complex phase does not occm' for smaller 
groups and, therefore, does not show up in the work of Macfarlane [100]. 



53 



principle. The determinant can be calculated by means of the formula^^ 

det U = Tq^ - tr - i (tr T^Y + tr + 8i To Pi T . (283) 

The sign of the last term is determined by our choice of the paravector space (specif- 
ically, our choice for the sign of the component cq = — 14). The explicit calculation 
of To in terms of the matrix L by means of eq. (283) turns out to be quite tedious. In 
view of the fact, that eq. (269) provides us with an explicit expression of Tq in terms 
of the matrix A we will not further pursue this subject here. The algebraic complex- 
ity of the calculation of Tq in terms of the matrix L seems to be related to the fact 
mentioned below from eq. (234) that the expression relating the antisymmetric ma- 
trix A to a sum of a finite number of powers of the matrix L is algebraically involved. 

The alert reader will have noticed that the method of Klotz [104] delivers twice 
as many results for the matrix U than the method of Fedorov and collaborators 
[235, 236] (Of course, the difference lies in the number of allowed values for the total 
phase of U only.) and than is required for the group SU{4) as the double cover of 
5*0(6; R). The explanation lies in the respective starting points of the two methods. 
In our version of the approach, the starting point of the method of Fedorov and col- 
laborators, eqs. (248), (249), is Z2-invariant for the elements of SU{A). On the other 
hand, eq. (276), the starting point of Klotz is invariant under the center Z4 of SU (4). 



5.3 From S0{6; M)/Z2 back to its double cover S0{6; M) 

We will now invert eq. (130), i.e., we will express the elements of the matrix L in 
terms of the elements of the (compound) matrix x = C2{L). This task occurs if 
one wants to find from the coefficients Xok^^ of ^ given Bogolyubov-Valatin transfor- 
mation the 50(6; M) matrix L generating it. We would like to discuss two possible 
lines of attack to this problem which, however, both have their deficiencies. 

The first method relies on eq. (277) [or eq. (280)] which can be inserted into 
eq. (248) to obtain the matrix L directly. As a first step, this method requires the 
(tedious) calculation of the complete matrix Xmn from the coefficients Xok^ [eq- 
(141) and its generalizations]. Furthermore, Tq is known from eq. (278) up to its 
phase only and its calculation entails the difficulties discussed at the end of Subsec. 
5.2.2. 

The second method to calculate the matrix L from the coefficients Xok^' takes a 
different route. Consider the following matrix elements (We also allow k,l = for 

32Cf. [72], Sec. 3.6, p. 42, eq. (3.62), [244], [98], p. 37 (reprint # 1: p. 55, reprint # 2: p. 948), 
eq. (44), [73], Sec. 65, p. 129, eq. (65.1), [245], Part II, §11, pp. 160/161, [206], Sec. 3, p. 819, eq. 
(40). Be aware of misprints and errors in the various versions of the equation as displayed in the 
references cited. 
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which Xok 01 vanishes individually.) 

XOk 01 = Z?/"" Xmk nl 

= {^iLpJ)Lki - {LI3j)^[l3iL)^ = LooLki - L^oLoi . (284) 

Of course, the choice of P = Pi [eq. (93)] involves a certain element of arbitrariness in 
expressing the matrix L in terms of x and is dictated here by physics considerations 
„„j „ — \^T„ „„ „, t — ™„ii,, T^Yite (assuming L[ -'- 

PjLPi 

XOk 01 + L^qLqi 



and convenience. We can now formally write (assuming Lqo 7^ 0) 

Xok 01 



'00 



(286) 



Note that Xo(-i) 0/ is given by 

~ Xo{-i) P1P2 ~ X(-i)o P1P2 
1 

~ ^ ^{-I)0kik2hl2 ^pip2giq2SiS2 

XOki mqi X0k2 mq2 X0l\ ns\ X0I2 ns2 (287) 

which is obtained by relying on the eqs. (128), (105) and = (i(_i)o = idid2d^dA,- 
The other components of xok oi can be read off from eq. (102). We now have to 
express L^o, Lqi in terms of xok oi yet. Let us start with Lqo- Taking the determinant 
on both sides of eq. (285) we find 

detL = (288) 

-^00 

where Xoo denotes the 5x5 matrix with the matrix elements Xok oi {k, 1^0). Taking 
into account the orthogonality of the matrix L (det L = 1) we finally obtain the 
relation 

iLool = (det xoo)^ • (289) 

The non-diagonal elements L^o, Lqi can be found by means of the following 
consideration. Taking into account the orthogonality relation LL^ = Ig "we find the 
following equation. 

Xofc°™ Xoi om = {Lko — Lqq Sko) {LiQ — Loo ^/o) + -^00 i^ki — Sko^io) (290) 
This equation can be transformed to read 

{Lko — Lqq 6ko) {LiQ — Loo ^/o) = Xofc°™ Xoi om — -^oo i^ki — ^ko^io) ■ (291) 
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Setting = / (7^ 0) we immediately find 

\Lko\ = ^Jxok"" Xok on.-Ll, = \/xofc°"^ Xok - (det xoo)^ • (292) 
Considering X^"^ok Xom 01 we obtain in an analogous way 

\Lok\ = V^X%fc Xom ok-Ll, = V^X^fc Xom ofc-(det Xoo)^ • (293) 

As an aside, we also find the following expressions for k ^ I 0) 

LkoLio = Xok"" Xoi om , (294) 

LokLoi = x°"^ok Xom. 01 , (295) 

which entail the following relations for the matrix xoo (no summation over k, I on 
the r.h.s.) 

n - 

Xok Xoi Om - Ski (det xoo) ' 

= y (xofc°™ Xok Om - (det xoo)^) (xo«°"' Xoi om - (det xoo)^) , (296) 

n - 

X ™ofe Xom 01 - Ski (det xoo) ' 

= y Xom ofc - (det xoo)^ (x% Xo™ oz - (det xoo)^) • (297) 

Collecting all results and inserting them into eq. (286), our final expression for 
the matrix L in terms of the matrix x (more specifically, in terms of its submatrix 
Xoo) reads 



Xofc 01 



\Lki\ = (det Xoo)' 

± y (xofc""" Xok Om - (det xoo)^) (x% Xom 01 - (det Xoo)^) 

(298) 

Already from eq. (130) it is clear that the matrix elements of the matrix L can be 
determined from the matrix x up to a minus sign only. The ±-sign in (298) reflects 
an uncertainty and, therefore, weakness of the applied method to reconstruct the 
matrix L from its second compound matrix x = C2{L). 
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6 Discussion 



In this section we collect the discussion of various aspects of the study performed 
in the preceeding parts of the article. We have postponed these considerations until 
now in order to streamline the presentation in the other sections. 



6.1 Bogolyubov-Valatin transformations 

To prepare ourselves for the following discussion we first have a look at the problem 
how the 15-dimensional irreducible representation of the Bogolyubov-Valatin group 
S'0(6;M)/Z2 operating in the biparavector space A^(^) decomposes into repre- 
sentations of smaller dimension under the subgroup chain 5*0(6; M) D 0(5; M) D 
0(4; M) D 0(3; M) D 0(2; R). This decomposition can be described by the graphi- 
cal representation given in Figure 2. 

^0(6; M) 15 



0(5; M) 



0(4; M) 




0(3; M) 3©3©3©1©1©1©3 




0(2; M) 2©1©2©1©2©1©1©1©1©1©2 



Figure 2: Decomposition of the 15-dimensional biparavector space /\^{Vq) under 
the subgroup chain of 5*0(6; M). 
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6.1.1 Linear Bogolyubov-Valatin transformations 

Now we will examine the well-known class of linear [0{2n = 4; M) C 5*0(6; M)]^^ 
Bogolyubov-Valatin transformations (Then, the only nonvanishing coefficients in the 



eq. (39) are A 



(1|0) ^{2|0) ^{0|2) 



A 



A 



and see how the decomposition of the bipara- 



vector space A^(^) structured under this class of canonical transformations. The 
6x6 matrix L G 5*0(6; M) then has the form (somewhat symbolically written) 




where the 4x4 matrix A G 0(4; ! 
(48), (49), and (A.30)]: 



The matrix A 



A 



/ Re 




Re 




Re 




Re 




\ 


Im 




Im 




Im 




Im 






Re 


^(l|0) 


Re 


^.(0|1) 
2 


Re 


Am 

2 


Re 


J0\2) 
2 




y Im 




Im 


^(0|l) 
"'2 


Im 


^(2|0) 


Im 


^(0|2) 
"'2 


/ 



(299) 

A ({A}) reads [cf. the eqs. 

- X^'l . (300) 



Here, x'^' denotes the matrix with row number k and column number m matrix 



elements Xfc' Eq. (299) can be found by relying on Subsec. 5.3 and on eq. (127). 

Consider now the linear space W (of operators in Fock space) which is the direct 
sum of the space spanned by the identity operator in spin space C4 (by the way, 
an invariant subspace) and the space of biparavectors A^(^6)- K is spanned by the 
basis^'^ 



4, 



za?^ A A a?', -^a? A a™ A a?,za? A af^ A a?', -za™ A af^ A a™, 



—a 



[1] 
1 ' 



■1" 



*2 ' 



-a 



[2] 
1 ' 

4" 



2 ' 
it'l - 

^-2 1 



.1" 



''2 1 



-a?Aa^\-a?Aa^\-a^^Aa^^^ 



(301) 



(I4, —ici A C2 A C3 A C4, 
ic2 A C3 A C4, —ici A C3 A £4, ici A C2 A C4, —ici A C2 A C3, 



■^■^Cf. e.g., [1], Sec. 2.2, p. 38. For further references concerning linear Bogolyubov-Valatin 
transformations see [49], p. 10247, below from eq. (16). 

■^''To achieve a more compact graphical display, we give here the Hermitian conjugate . Be 
aware of the minus signs and note that in the last line for simplicity we have used a somewhat 
imprecise notation [cf. eq. (129)]. 



58 



— Ci, — C2, — C3, — C4, 



-Ci A C2, -Ci A C3, -Ci A £4, -C2 A C3, -C2 A q, -C3 A £4, ) 

Xt it it it ;t 

U, C(_1)0, C(_;^)^, C(_;^)2' C(_l)3, C(_;^)4, 



jt xt ;t xt xt ;t xt ;t xt xt 
^015 ^025 ^03' "^045 ^125 ^13' ^14' ^23? C24, C34 



U4,<^ I • (302) 



Also defining'^^ 

6t = (l4,-z6flA6flA6WAS?l, 



zfefl A foW A -S? A feW A fel^zfeW A fof A fct^ -zfo™ A fef^ A 6^, 

-fotil _5[2] 
f^i ) "2 ' "1 ' "2 ' 

-6WA6f],-6WA6W,-6WA6[^ 

-6f]A6W,-6riA6[^-6|lA6[^]) (303) 
U4, — ic/i A (i2 A (is A (i4, 

id2 A c/3 A (^4, — iiii A ^3 A ^4, idi A ^2 A ^4, — idi A ^2 A ^3, 
—di A (i2, — A d-s, —di A ^4, —d2 A 1^3, — ^2 A d^, —d^ A (i4, 

J t X t X t J t X t 

I4, d (_i)0; (-l)l' (-1)2' ^ (-1)3' (-1)4' 

it X t X t X t X t X t X t X t X t X t 
d Qi^d Q2id 035 d Q4, d d 13, 14, d 23, 24' 34 

li,D 1 (304) 



^^Below we have not indicated that the wedge product relates to a different vector space than 

b d 

in eq. (301), i.e.: A = A = A. 
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we can write any two-mode Bogolyubov-Valatin transformation as a basis transfor- 
mation in the linear space W. 

b = Ai{X}) a (305) 

For a linear Bogolyubov-Valatin transformation, the matrix A assumes a block di- 
agonal form. It reads [cf. eq. (127)]: 

^({A}) = dis,g[CoiA),C,{A),Cs{A^)\C,iA).C2iA)] (306) 

Here, Ci{A) denotes compound matrices of the matrix A (cf. Appendix C). Let us 
denote the subspace in which the compound matrix Ci{A) operates by Wi (W = 
Wq © W4 © W3 © Wi © 1^2)- With the usual relations for compound matrices and 
the orthogonality condition for the matrix A we can write eq. (306) as: 

^({A}) = dmg[l,detA,{detA)A,A,C2{A)] (307) 

For the indices of the compound matrices we have applied here the usual convention 
of the lexicographical order of the indices. Note that the above representation of A 
in terms of compound matrices requires (and, therefore, induces) a certain order in 
the sequence of the basis elements of the linear subspaces W2, W3 invariant under 
the 0(4; M) transformation [cf. eqs. (301), (303)] once the sequence in the subspace 
Wi (spanned by d^^"^' = Ci, df^ = £2, 0.2' = C3, = C4) has been chosen. 



6.1.2 Nonlinear Bogolyubov-Valatin transformations for one fermion 
mode 

To make further contact with previously known results, we consider now the case of 
a full nonlinear 50(3; M) Bogolyubov-Valatin transformation of one fermion mode 
(say, with mode number 1) considered in 



(308) 

(309) 
(310) 



k = Bi ({A}; {a}) = Af ^ a, + aS^I°) a\ + A^^l^^ (^a\d^ - ^ I4) 
The coefficients Xf^^\ aS^'°\ A^^'^^ obey the equations (cf. [49], p. 10247) 



|A1 



(1|0)|2 



{0|1)|2 



2 ' ^ 



(1|1)|2 



4 AS^I°)Af + (aS^W 



. 



For this case, the 6x6 matrix L G 5*0(6; M) has the form [cf. eq. (127)] 



P(-1)(-1) 











^(-1)3 


^(-1)4 





Qoo 


Qoi 


Q02 











Qio 


Qii 


Q12 











Q20 


Q21 


Q22 








^3(-l) 











-P33 


P34 


^4(-l) 











P43 


P44 
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(311) 



where the 3x3 matrices P, Q ^ 5*0(3; M). The orthogonal matrix Q is related to 
the orthogonal matrix A ({A}) used in eq. (23), p. 10248 in [49] the following way 



Am) 



Re 
Im 



Re 
Im 



Re 
Im 



^ Im(^fi:(o|i)fi:(i|i)) Im ( ) Im ( 

C2 (g) . 



Here, 



K 



(1|0) 



ii) 



A(oii) + x'^m ^ 



Relying on the equation (for the notation applied cf. Appendix C) 



Ci {QY 



(312) 

(313) 
(314) 
(315) 

(316) 



where 




(317) 



we quickly find 

Q 



eAi{\}) e 
( Im 

-lmK(ili) 
\ Re 



-Im 

Im 
-Re 



-Im /^(ilo) 
Re / 



(318) 



The orthogonal 3x3 matrix P must be a representation of the special orthogo- 
nal transformation given by the matrix Q (in view of the group properties of gen- 
eral Bogolyubov-Valatin transformations). The simplest choice consists in setting 
P = I3, however, also putting P = Q (or setting it equal to some equivalent repre- 
sentation matrix) is possible. At this point we are mainly interested in the question 
of whether a full nonlinear Bogolyubov-Valatin transformation of the fermion mode 
with mode number 1 necessarily induces a related (nonidentical) transformation of 
the fermion mode with mode number 2 in such a way that the combined transfor- 
mation is canonical. One can now convince oneself of the fact that no admissible 
choice of the matrix P leaves the fermion mode with mode number 2 unchanged. 
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For the "minimal" choice P = I3 we find 
k = B2({A};{a}) 



ai 



(319) 



Of course, for a hnear Bogolyubov-Valatin transformation of mode 1 (i.e., A^^'^"* = 0) 
the two fermion modes decouple. 

This consideration demonstrates the truly nonlinear character of general Bogo- 
lyubov-Valatin transformations. A nonlinear Bogolyubov-Valatin transformation 
for one mode necessarily requires us to perform corresponding changes for the other 
mode(s) in order to maintain the correct CAR. The picture of independent quasi- 
particles characteristic for linear Bogolyubov-Valatin transformations no longer ap- 
plies in general. 



6.1.3 Nonlinear Bogolyubov-Valatin transformations 
on the ansatz of Fukutome and collaborators 



Some comment 



One of the earlier studies of nonlinear Bogolyubov-Valatin transformations is con- 
tained in an article be Fukutome and collaborators [31]. Our study of these transfor- 
mations allows us to see the ansatz used in ref. [31] from a different perspective. For 
the interested reader (who, for easy reference, is advised to have a copy of the article 
[31] close at hand), we would like to mention that our eq. (105) bears some resem- 
blance to eq. (2.11) on p. 1557 in [31] specialized to the case of two fermion modes. 
However, while in part it corresponds to our more general equation (105) the eq. 
(2.11) in [31] does not adequately reflect in its display the general group theoretical 
structure of nonlinear Bogolyubov-Valatin transformations. To see this note that 



the case considered in ref. [31] corresponds to the choice xlf^* = 0, Xk^* = in our eq. 
(102). Consequently, under the reduction 5*0(5; M) C 5*0(6; M) considered in [31] 
the 15- dimensional irreducible representation of the nonlinear Bogolyubov-Valatin 
transformation [equivalent to 5*0(6; M)/Z2] operating in biparavector space A^(^) 
splits into two irreducible representations of 5*0(5; M) - a 5-dimensional representa- 
tion [operating in the space of biparavectors with indices k = (—1) and / ^ (— 1)] 
and a 10-dimensional representation [operating in the space of biparavectors with 
indices k,l 7^ (—1) ]; cf. Figure 2. The 5-dimensional linear space (of operators in 
Fock space) is spanned by products of three or four different fermion creation and 
annihilation operators while the 10-dimensional space (of operators in Fock space) 
is spanned by the fermion creation and annihilation operators themselves and all 
possible products of two (different) of them. The unit operator in Fock space does 
not transform under nonlinear Bogolyubov-Valatin transformations and it spans an 
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invariant subspace of the 16-dimensional space of linear operators in Fock space. 
Consequently, eq. (2.11) on p. 1557 in [31] must be the above characterized 5- 
dimensional representation in disguise and not a transformation in paravector space 
[associated with the Clifford algebra C(0,4)] as its display suggests. 

6.1.4 The group of nonlinear Bogolyubov-Valatin transformations as the 
adjoint representation of the group SU (2" = 4) 

Recalling eq. (127) 

USmU^ = Xa/" (320) 

we recognize ([275], Vol. II, Chap. 11, Sec. 5, Theorem III, p. 417) that the matrix 
(i.e., the transposed of x) is the representation matrix of the group element U of 
the Lie group S'f/(2" = 4) taken in the adjoint representation, i.e., we can write 

= [X{^)f = Ad(U) . (321) 
Consequently, we find [cf. eq. (127)] 

Ad(U) = C2{L^) (322) 
where the matrix L = L (U) G 5*0(6; M) is given by [cf. eq. (249)] 

Lki = I tr (f , U) . (323) 

Clearly, the adjoint representation of the Lie group SU{2'^ = 4) is not faithful 
because S'[/(2" = 4) posesses a non-trivial center (i.e., Z4). Let us finally note that 
the adjoint representation of the group SU{3) has been discussed in [246]. 

6.2 Biparavectors and parametrizations of SU{4) and 
50(6; M) transformations 

Low-dimensional unitary and orthogonal matrices play an essential role in many 
investigations in theoretical physics and beyond. The choice of representation for 
them can significantly facilitate or inhibit their use. In the mathematical and physi- 
cal literature many different parametrizations can be found. Concerning the unitary 
matrix U G SU (4) we are of course interested in representations in terms of the bi- 
paravectors c M, i-e., in terms of fermion creation and annihilation operators (For a 
related discussion see [247, 248].). We have discussed in Sec. 5 such a representation 
of SU{4:) matrices in terms of biparavectors [cf. eq. (226)]. In this subsection, we 
would like to supplement this discussion by some further comments. 



63 



6.2.1 The representation of S'f/(4) transformations by Ostlund and Mele 

An interesting representation of tlie unitary operator U related to eq. (226) (i.e., in 
terms of biparavectors) has been considered by Ostlund and Mele [43]. Define the 
following states (We have changed somewhat the numbering in comparison with ref. 
[43].): 

^(-1) = |2)(2,1) 

*o = |0) , 
*i = |l)i = 

*2 = |1)2 = 

Then, one defines the elements of the matrix rh (= m, in the notation used in [43]) 
by means of the relation 

^ij = . (328) 

It is then useful to recall the relation 

|0)(0| = ^o^S = (I4 - aiai) (I4 - 4^2) (329) 

to explicitly construct the matrix m = m in terms of the creation and annihilation 

operators d|., (cf. eq. (1) of ref. [43]). Consequently, the matrix elements rhij of 

the matrix m are linear combinations of the operators Cm (and transform under 
S'0(6;M)/Z2 transformations correspondingly). For the explicit relations see Ap- 
pendix F. 

One can now write {Xij G C) 

2 

U = U(X) = J2 ■ (330) 

A remarkable property of the representation of the unitary operator U in terms of 

the matrix rh is represented by the following equation where XY denotes the matrix 
multiplication of the two 4x4 matrices X and Y. 

U(X)U(r) = U(XF) (331) 

Taking into account the relation [U(X)]^ = U(X'I") (eq. (3) of [43]), the unitarity 
condition for the unitary operator U is then found to read [This is another version 
of the eqs. (227), (228).] 

XX^ = I4 . (332) 

Consequently, eq. (330) provides us with a representation of the unitary matrix U 
by means of some, in general other, unitary matrix X. 



ala\\0) 



al\0) , 
4|0) . 



(324) 
(325) 
(326) 
(327) 
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6.2.2 Exponential representation 

Next we want to have a short look at some other version of the parametrizations 
of the unitary [U G SU{4)] (and orthogonal [L G 5*0(6; M)]) transformations we 
have applied. Besides the representation we are using [see eq. (226)] of most interest 
appears to be an exponential representation of the unitary matrix U 

U = To I4 + T'' Im = To I4 + ^ T^'"^' L,^, 

= e^^{y''lM) = exp Q L,^,^ (333) 

in terms of a real antisymmetric matrix V (remember: Cj^/ = —c m)- Formally, 

V'Im = \ V^^^^ lm,m2 = InU . (334) 

There is a comprehensive mathematical literature on the subject of the logarithm 
of a matrix and its subtleties. Restricting our attention to the main sheet of the 
logarithmic function, a calculation of the matrix V in terms of Tq and the matrix T 
could rely on the approach presented in [249], in particular Sec. 2.C.3, p. 021108-5 
(also note [250]). However, we will not pursue this calculation here. 

Concerning some of the subtleties of the matrix logarithm it is useful to recall 
that any n x n unitary matrix operating in a complex space of n dimensions can be 
represented as a 2n x 2n orthogonal matrix operating in a real space of dimension 2n. 
Consequently, the problem of finding the exponential representation of the unitary 
matrix U is closely related to the problem of the exponential representation of the 
orthogonal matrix L in terms of an antisymmetric matrix B 

L = . (335) 

For our case of the 6x6 matrix L, by virtue of the Cayley-Hamilton theorem the 
antisymmetric (6 x 6) matrix B must be a linear combination of the matrices A^, 
k = 1,3,5, in the vicinity of the unity transformation L = Ug- For the subject of 
the logarithm of orthogonal matrices we refer the reader to [251-253]. Finally, we 
want to mention that for bi (para) vectors two different notions of the exponential 
function exist (see [254], [56], Sec. 17.3, p. 221, both editions). But, we will not 
further discuss this subject here. 



6.2.3 On Cayley-Klein parameters for 5*0(6; R) transformations 

In the light of the equations (246), (247) (Subsec. 5.1) reflecting the relation between 
50(6; M) and SU (4) transformations it seems to be interesting to add a comment on 
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the introduction of generalized Cayley-Klein parameters for the orthogonal trans- 
formations L in the Euclidean space Mg- For ^ good introduction to Cayley-Klein 
parameters in M3 we refer the reader to [255], §4.5, 1. ed. p. 109, 2. ed. p. 148 (Note, 
that in the third edition [256] this material is no longer included.). Leaving aside the 
problem of the signature of the (pseudo-) Euclidean space on which the (pseudo-) 
orthogonal transformations operate, Cayley-Klein parameters for (pseudo—) orthog- 
onal transformations in a 4-dimensional space have been studied by Cayley already 
[257]. They have later been discussed in [258-262], [263], also see [264], Mathemati- 
sche Zusatze und Erganzungen, Sec. 17, pp. 806-813, [89, 265-269]. Corresponding 
studies of Cayley-Klein parameters for (pseudo-) orthogonal transformations in five 
dimensions can be found in [270-272]. A five- dimensional space seems to be the 
space of largest dimension to which the standard Cayley-Klein parameters can be 
generalized. From our perspective, the reason is that for orthogonal transforma- 
tions L G 5*0(5; M) the 15-dimensional biparavector space A^(^6) decomposes into 
two invariant subspaces. One subspace is a 5-dimensional space in which the five- 
dimensional transformation L acts linearly. The second invariant subspace is a 
10-dimensional space in which the 5-dimensional transformation L acts via its sec- 
ond compound matrix exactly in the same manner as an orthogonal transformation 
L G 50(6; M) is acting in the full 15-dimensional biparavector space A^(^6) [cf- 
eq. (127) and Fig. 2]. Consequently, generalized Cayley-Klein parameters for or- 
thogonal transformations in Mg must be based on the spin bivector space 
according to the eqs. (246), (247), rather than on eq. (127) related to biparavector 
space A^(^6)- 



6.3 Biparavectors as bi(para)vectors of the spin bivector 



Here, we would like to extend somewhat the discussion performed in paragraph IV 
of Subsec. 3.2.3. By means of the eqs. (242), (245) (Subsec. 5.1) one finds for the 
decomposable/simple biparavectors [cf. eq. (23); mi 7^ rri'if'^ 




c. 







(336) 



For the related references see footnote 23 to eq. (242) in Subsec. 5.1. 
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In view of eq. (336), we can write 

/\\V,) = /\' (Ai(C4)) = , (337) 

but clearly Vg 7^ Mg (The vectors of the spin bivector space Mg are antisymmetric ma- 
trices while one direction in the paravector space Vg is given by the unit matrix I4.). 

Eq. (336) tells us that the decomposable/simple biparavectors cm of the paravector 
space Vg can also be understood as decomposable/simple bi (para) vectors of the spin 
bivector space Mg ^"^ . In other and more plain words, any (4 x 4) Dirac matrix and 
any of their (multiple) products can be given as the product of two antisymmetric 
(4 X 4) matrices that exhibit a certain generalized anticommutation relation. This 

insight into the structure of decomposable/simple biparavectors cm has in part been 
spelled out by Haantjes, [144], p. 51, stelling 5 [proposition 5], and a corresponding 
comment can also be found in ref. [94], §14, p. 269, above of eq. (14.8). While we 
have concentrated our study on the structure of S'0(6;]R) transformations in bipa- 
ravector space A^(^) Haantjes [144] and Buchdahl [94] reflect the decomposability 
property from out the construction of the paravector s/biparavectors in terms of the 
product of spin bivectors [cf. eq. (242)], i.e., their representation in terms of prod- 
ucts of antisymmetric matrices. For the special choice discussed in Appendix D this 
product representation of the biparavectors corresponds to a product representation 
in terms of two commuting sets of quaternions H. Another product representation 
(in terms of matrices called Sa and Da) has been discussed by Eddington [68], [72], 
Chap. Ill, Sees. 3.1, 3.2, pp. 34-36. Nikol'skii [273] has shown that this product rep- 
resentation is related to transformations of tetrahedral coordinates (For a related 
discussion concerning the Dirac equation also see [274].). 



6.4 Miscellaneous 

Finally, we would like to add some comments related to the choice of the metric (in 
particular its signature) in spin bivector space. According to the eqs. (246), (247) 
(Subsec. 5.1), unitary transformations [g SU{4:)] in the spin space C4 correspond to 
orthogonal transformations [g S0{Q; M)] in the spin bivector space Mg. On the other 
hand, it has been known for a long time that quadriquaternions (or, equivalently, the 

biparavectors cm) are related to conformal transformations [83]. To arrive at these 
one must change from the spin bivector space Mg to the related spin bivector space 
Ms^i. This change in metric can be achieved, for example, by means of the transition 

— y i r_]^ [entailing, e.g., (r_^] = — 14]. The 5*0(5, 1;M) transformations in 

^''We tend to use the term 'biparavectors' of spin bivector space here because their definition 
involves a Hcrmitian conjugation in one of the two factors in analogy to the biparavectors defined 
in eq. (23). 
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1^5,1 then stand in correspondence to SU*{4:) transformations in spin space^^. It 
is then possible to represent the pseudo-orthogonal 50(5, 1;R) transformations in 
Ms i by means of conformal transformations in a subspace M4 C Ms^i. However, 
the corresponding transformations in spin space then are no longer unitary and 
do not correspond to canonical transformations in which we are interested. From 
a particle physics perspective, transformations not in Euclidean space M4 but in 
Minkowski space Ma^i are interesting. The Minkowski space is then to be embedded 
in a six-dimensional pseudo-Euclidean space and pseudo-orthogonal 5*0(4, 2; M) 
transformations in M4 2 then correspond to pseudo-unitary SU{2, 2) transformations 
in spin space. We will not further dwell on this subject here but confine ourselves to 
some comments on the related literature: For a good, general mathematical account 
of conformal transformations see [277], Chap. 11, p. 480. Ref. [278] covers the subject 
from a physical perspective (including many references, also see [279]). For articles 
dealing with various aspects of the above relationship that are interesting within 
the present context see [95, 107, 142, 143, 280, 281]. This subject has also been dis- 
cussed, in part recently, from a twistor perspective (see, e.g., [104, 282, 283]). Leav- 
ing aside special conformal transformations and dilatations in the pseudo-Euclidean 
space 1^3,1, the relevance of the restricted class of 5*0(4, 1;M) transformations in 
R42 for massive (4-component) fermions (Dirac particles) has been discussed in 
the physics literature (Dirac equation in 4 -|- 1 dimensions) [270, 284-287]. Also, 
(generalized) Foldy-Wouthuysen transformations in spin space fit into the latter 
framework [288, 289]. A discussion of the algebra of Dirac matrices related to the 
group 50(3, 2; M) can be found in [290-292]. Finally, for the sake of completeness, 
we also want to mention the papers [293] and [294-298] which have attempted to 
address certain problems related to the general spin space transformations we have 
studied. 



7 Conclusions 

Having reached the end of the present study of nonlinear Bogolyubov-Valatin trans- 
formations for two fermion modes, it seems fair to say that even this apparently 
simple and basic case has revealed its complexity and broad significance. Among 
the more general aspects, we would like to point out that unlike linear Bogolyubov- 
Valatin transformations which can be performed for each mode of quasifermions in- 
dependently and separately, nonlinear Bogolyubov-Valatin transformations exhibit 
the truly nonlinear nature of the canonical fermion anticommutation relations be- 
cause any nonlinear Bogolyubov-Valatin transformation - even within one fermion 
mode - immediately impacts the other mode(s) and, in general, the (linear) picture 
of independent quasifermions no longer applies. 



^^For the notation see, e.g., [275], Vol. 2, p. 392, Table 10.1. S'[/*(4) is the nonconipact version 
of SU{A) related to it by the "unitary trick" of Weyl [276]. 
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From a mathematical point of view, nonlinear Bogolyubov-Valatin transforma- 
tions are accompanied by transformations of multivectors of different grades and, 
therefore, can be viewed as concrete realizations of supersymmetric transforma- 
tions. The preservation of the canonical anticommutation relations under nonlinear 
Bogolyubov-Valatin transformation is closely connected to the decomposability of 
the biparavectors generating the Clifford algebra C(0, 4). This is the Clifford algebra 
that is associated with the creation and annihilation operators of the two fermion 
modes. From a somewhat different perspective, the preservation of the canonical 
anticommutation relations under nonlinear Bogolyubov-Valatin transformation can 
be described in other words as being closely connected to the fact that Dirac ma- 
trices can always be written as the product of two antisymmetric matrices. This 
is a fact that was first recognized in the 1930's already by Dutch mathematicians 
(Schouten, Struik, Haantjes). 

Looking to future generalizations of the present work to more than two fermion 
modes it should be said that one should expect to meet mathematical structures 
which have not let their presence in the study of one and two fermion modes be 
known. The present investigation of nonlinear Bogolyubov-Valatin transformations 
for n = 2 fermion modes exhibits a number of features that are related to fact that 
the group SU{2'^ = 4) is the double cover of the group 5*0(6; M) [For one fermion 
mode SU(2) is related to S'0(3;M).]. For n > 2 fermion modes no such relation 
of the group SU (2") to any of the orthogonal groups exists. This fact necessarily 
will raise its head in any future study of nonlinear Bogolyubov-Valatin transforma- 
tion for more than two fermion modes. This will lead to a further increase in the 
complexity of the problem to be studied but certainly also to further interesting 
mathematical and physical insight into the nature of fermionic systems. 
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Appendix A 



We define the coefficients [cf. Sec. 2, eqs. (39), (40)] 
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This in turn implies: 



A 



(0|1,2) 



A 



(1,2|1) 



A 



(1|1,2) 



A 



(1,2|2) 



A 



(2|1,2) 



A 



(1,2|1,2) 



1 

2 

1 
2 

1 
2 

1 
2 

1 
2 

1 
2 

1 
2 

1 
2 



K 



(1|0) 



..(Oil) 



+ K 



(1,2|2) , ■ (2|1,2)\ 



+ IK, 



(2|0) , - 



(012) 



.(1>2|1) 



(1|1,2)A 



-K 



,.(2|1) (1|2) 



IK 



(1,2|0) 



{0|1,2) 



..(2|1) ^ fl|2) 



,.{1,2|0) 



IK 



,.{0|1,2) 



(2|1) 



(112) 



K 



{1,2|1) 



.(2|1) 



.(1|2) 



tK 



(1,2|0) 



(0|1,2) 



(1,2|0) 



IK 



(0|1,2) 



(1|1,2) 



— K 



(1,2|1) , .(111,2) 



— K 



(1,2|2) 



— IK 



(2|1,2) 



..(1,2|2) 



— IK 



(2|1,2) 



2 K 



(1,2|1,2) 



(A.17) 

(A.18) 

(A.19) 

(A.20) 

(A.21) 

(A.22) 

(A.23) 

(A.24) 
(A.25) 



(A.26) 
(A.27) 
(A.28) 
(A.29) 

Furthermore, the following relations apply [cf. eqs. (48), (49) and (55); k = 1,2] 
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Appendix B 

To work out the anticommutator (and commutator) of the operators dk (see Subsec. 
3.2.1) it turns out to be useful to calculate the following products first (On the r.h.s. 
terms with more than four factors in the wedge product have been omitted as they 
vanish in the case under consideration.). 

c.di = - xfL) l4 - xf 1,,.) £ ' + xF' A - i xfls,t) ^ A c 

1 [91 (r si . . . 1 



- ^ xP^ c, A a. A a. - ^ xlt,,,.) . a ^ A a * A a - 

+ ^ xF A a, A a, A Q (B.l) 

p ^ q., the r.h.s. is antisymmetric in p and q: 

" - 3 [2] n [1] - , [1] - [3] ■ - t 

- xf £p A a ^ + xf\p,s) A a ^ 



2! 



Xl \p,,,t,u) Ac + Xl ' ' CpACgACr 



- ^ xl ^ Cp A c ^ A c ' + - x; |p,,,t) ^ A c ^ A c ^ 

^ ^'''^ ap A a, A a, A a, . (b.2) 
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n 7^ p 7^ g, the r.h.s. is completely antisymmetric in n, p and q: 

— - i _ [3] • H _ [2] - I [2] _ [2] 

CnCpCqai — Xl (n,p,q) ^ ^4 X; (p,g) C„ + Xl (n,q) ^/ (n,p) 



+ Xf\n,q,t) iCpAc' - Xf\n,p,t) ^ ^q N 

- c„ A Cp A c ^ + xf (p,,) c„ A c, A c ^ 

- Cp A A c ^ + xf^ c„ A Cp A A . (B.3) 

m ^ n ^ p ^ the r.h.s. is completely antisymmetric in m, n, p and q: 

- - - - i _ [4] • -n _|_ [3] • - _ [3] • - I [3] 

CmCnCpCqai — X/ (m,n,p,g) ^ -"-4 + (n,p,g) ^ '^m X/ (m.,p,q) ^^n + Xl {m,n,q) ^ '^P 

[3] . [2] ^ -X [2] ^ 

~ Xl {m,n,p) ^ '^q ~ Xl (p^q) Cm /\ Cn + Xl (n,q) ^ 

[2] ^ ^ [2] ^ ^ [2] ^ ^ 

~ Xl (n,p) Cmf\Cq — Xl ^ + X/ "^n ''^ "^9 

~ {m,n) Cp ACq — x\ \q) Cm A Cn A Cp + xl \p) A C„ A Cg 

- xl'(„) CniACpACq + xl'(™,) c„ A Cp A Cg . (B.4) 

In deriving the above relations we have used the Clifford algebra relations (x, are 
vectors and A is an arbitrary multivector.). 

X A = xJA + X A A (B.5) 

and (y/c denotes a vector to be omitted.) 

p 

xj (yi A... Ayp) = (x-y,,)(yi A... AyfcA... Ayp) (B.6) 

k=l 

(Cf., e.g., ref. [f52], Chap. 1, p. 9, eq. (1.33). Note, that Hestenes and Sobczylc use 
a somewhat different notation for the operation of left contraction J.). 



Appendix C 

To make the present paper self-contained we repeat in this Appendix some defini- 
tions and mathematical results concerning compound matrices given in [299], Subsec. 
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II. A, p. 5420 and Appendix A, pp. 5443/5444. 



We introduce x matrices A^^'^) (/c = 1, . . . , n) by writing (choose li < I2 < 
. . . < Ik, nil < m2 < . . . < ruk) 

.(2k) _ .{2k) .^.x 

^LM ~ ^h...lk,mi...mk \^-^) 

(We identify the indices L, M with the ordered strings li . . . 1^, mi . . . rrik-) or, more 
generally (not requesting li < I2 < . . . < h, nii < 1122 < ■ ■ ■ < mi) 

^LM = sgn [aaih, Ik)] sgn [(Tfe(mi, . . . , mfc)] A[3^,mi...mfc • (C-2) 

The indices L, M label the equivalence classes of all permutations of the indices 
h, . . . ,lk and mi, ... , m^, respectively, and cxa, db are the permutations which bring 
the indices li,mi {i = l,...,k) into order with respect to the < relation (i.e., 
c^aik) < CTaih) < ... < cTaik), o^b("2i) < ^^(ms) < ... < ab{mk)). The matrix 
elements of the matrix A^^*^^ are arranged according to the lexicographical order of 
the row and column indices L, M (We identify the indices L, M with the ordered 
strings aa{li) ■ ■ ■ (Taih), crbirrii) . . . ab{mk), respectively.). 

We also define a set of (dual) (^) x matrices A^^'^)* {k = 1, . . . , n) by writing 

^i2k)* ^ g(k) i^{2k)T g{k)T j^Q3^ 

where the (^) x matrix S^^'^ is defined by 
consequently, 

g(k)T ^ f^_^Y"--k)k g{n-k) (C.5) 

(Quite generally, for any (^) x (^) matrix B we define B* by B* = £^'''>B'^£^^^^ .). It 
holds (Ir is the r x r unit matrix) 

gik) g{k)T ^ ^^^^ ^ ^C_Q) 
g{k)T g(k) ^ ^^^^^ ^C.7) 

We can now give some formulas for compound matrices'^^. Let B, D be n x n 
matrices. The compound matrix Cfc (B), < k < n, is a (^) x (^) matrix of all 
order k minors of the matrix B. The indices of the compound matrix entries are 
given by ordered strings of length k. These strings are composed from the row 
and column indices of the matrix elements of the matrix B the given minor of the 



■^^In the context of projective geometry, these matrices and their elements often are referred to 
as Pliickcr-Grassmann coordinates. 
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matrix B is composed of. Typically, the entries of a compound matrix are ordered 
lexicographically with respect to the compound matrix indices (We also apply this 
convention.). The supplementary (or adjugate) compound matrix C "^'^ (B) (some- 
times also referred to as the matrix of the kth cof actors) of the matrix B is defined 
by the equation (cf. eq. (C.3)) 

C^-'{B) = C„„fc(Br • (C.8) 

The components of the supplementary compound matrix C^~^ (B) can also be de- 
fined by means of the following formula (here, li < I2 < ■ ■ ■ < Ik, ^1 < ^2 < ■ ■ ■ < 
mk] [300], Chap. IV, §89, p. 75, [301], Chap. 3, p. 18)^°. 

C {B),M = ■ ■ ■ <^ det B (C.9) 

This comparatively little known definition of (matrices of) cofactors (supplemen- 
tary compound matrices) is essentially due to Jacobi [302], §10, p. 301, p. 273 of 
the 'Gesammelte Werke', p. 25 of the German transl. (also see the corresponding 
comment by Muir in [303], Part I, Chap. IX, pp. 253-272, in particular pp. 262/263). 

For compound matrices holds (1^ is the r x r unit matrix, a some constant) 

Ckialn) = a' . (C.IO) 
Important relations are given by the Binet-Cauchy formula 

Cfc(B)Cfc(D) = Ck{BD) (C.ll) 
from which immediately follows 

Cfc(B-i) = Cfc(B)-' , (C.12) 

the Laplace expansion 

Cfc(B)C"-*^(B) = C"-*^(B)Cfc(B) = 

Ck{B)Cn-k{BY = Cn-k{BTCk{B) = detBl(„), (C.13) 
Jacobi's theorem (a consequence of the eqs. (C.13) and (C.12)) 

Ck(B-') = -^C'-'^iB) = -^Cn-kiBY , (C.14) 
^ ^ det B ^ ^ det B ^ ^ ' ^ ^ 

and the Sylvester- Franke theorem 

detCfc(B) = (detB)(^-i) . (C.15) 



4"Notc, that in the eqs. (A.2), p. 5443, (31)-(33), p. 5421 in [299] the indices I and m on the 
r.h.s. should be interchanged to correct the display of these equations. 
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Compound matrices are treated in a number of references. A comprehensive 
discussion of compound matrices can be found in [304], Chap. V, pp. 63-87, [305], 
Chap. V, pp. 90-110, and, in a modern treatment, in [306], Chap. 6, pp. 142-155 
of the Enghsh translation. More algebraically oriented modern treatments can be 
found in [147], Part I, Chap. 2, Sect. 2.4, pp. 116-159, Part II, Chap. 4, pp. 1-164 
(very thorough), [307], Chap. 7, Sect. 7.2, pp. 411-420, and [308], Vol. 3, Chap. 2, 
Sect. 2.4, pp. 58-68. Concise reviews of the properties of compound matrices are 
given in [309, 310]. Also note [311]. 



Appendix D 



In this Appendix, we first give a concrete example for the complex antisymmetric 
matrices introduced in Subsec. 5.1: 
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They represent spin space bivectors (g A^(^4))- For k = 0,1,2, the relation 

+ 

= apphes, consequently, these three matrices span the eigenspace of the 

Hodge operator to the eigenvalue -1 [cf. eq. (236)]. For k = —1,3,4, the relation 

+ 

= — applies and, consequently, these three matrices span the eigenspace of 

the Hodge operator to the eigenvalue 1 (For a related discussion see [312].). The 

+ 

antisymmetric matrices with indices /c = 0, 1,2 and k = —1,3,4 are related to 
a matrix representation of two commuting copies, I'^, J'^, K'^ (I'lJ'i = K'^, etc.) and 
I2, J2, (I2J2 = 1^2, etc.), of the system of quaternions ([262], p. 329, eq. (5)). 
We attach here a prime to the quaternionic units in order to distinguish them from 
the different choice made in the eqs. (29)-(34) in Subsec. 2.2. With some hindsight 
we set 



I'l = 


i F2 = 


i ^2 > 


(D.7) 


J'l = 




-i F^ , 


(D.8) 


k; = 


+ 

i Fq 


^ To , 


(D.9) 


I'2 = 


+ 




(D.IO) 


j'2 = 


+ 

r4 = 




(D.ll) 


K'2 = 




- ^(-1) • 


(D.12) 



Note, that these relations are valid only if the matrices F^ with the indices = 0, 1, 2 
and k = —1,3,4, respectively, span the two eigenspaces of the Hodge operator to 
the eigenvalues —1 and 1. This will not be the case in general. Taking into account 
the above identifications we find the following relations of the quaternionic units to 

- coi = ci (D.13) 

- do2 = C2 (D.14) 

- C12 = C1C2 (D.15) 
C(_i)4 = i C1C2C3 (D.16) 

- C(„i)3 = i C1C2C4 (D-17) 

related approach has been used in [180], See. 4.1, pp. 99-102, however, without relating the 
quaternionic units to the bivectors of spin space. The approach in [180] has rather some similarities 
with the one used by Milner [81]. Incidentally, a paper by Rau [176], p. 4, eq. (14), contains some 
related considerations. 



the elements of the biparavector space 

+ 



I'l 


= J'lK'i 


= -ToFi 
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= k;i; 


= -foF2 




= I'lJ'l 


= - ^1^2 
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= J2K'2 
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= K'2r2 
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— _ r 
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K = I2J2 = -rof2 = -du = -c-ih (D.18) 

For the generators of the Chfford algebra C(0,4) we obtain^^ 



and, consequently, 



Cl = 


Cl 


= I'l, 


(D.19) 


C2 = 


C2 


= Ji , 


(D.20) 


C3 = 


C3 




(D.21) 


C4 = 


C4 


= ^ , 


(D.22) 


C(-l) 






(D.23) 



The above identifications are up to minor detail the same as given by Schouten [75], 
p. 106, eqs. (5), (6), and Lemaitre [79], p. 170 (For related discussions also see [72], 
Sec. 2.3, p. 23 and Sec. 3.8, p. 47, [73], Chap. VI, § 53, p. 108, [81].). 



Appendix E 

Here we calculate traces of products of the operators [For related results see 
[100], p. 136, eqs. (2.12)-(2.21).]. These expressions are necessary for the discussion 
in Subsec. 5.1 [see eqs. (248)-(250)]. We find 



- tic I = trl^FoF, j = 4 5oi (E.l) 

which is a special case of 

tT{l\l^ = tr(f,F,) = A 6,1. (E.2) 

Furthermore, we have 

- tr^C/C^C„^ = tr (^FqF^ FmF„ j = 4: {601 6mn - ^Om ^In + ^On ^Im) (E.3) 

which is a special case of 

tl (^C,C I C^Cr^ = tl (Vf^T I TmTn^ = 4: (Ski 5„m - Skm Sin + 5kn 5lm) , (E.4) 



"^^This choice corresponds to the Cartan extension of the Chfford algebra C(0, 2) generated by 
61, 62, cf. [313], Sec. 6.3, p. 80, eq. (6.26). 
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and 



tr I CiC^CnCpCq J — tr ( Fq VmXn^p^q 



4 ^Olmnpq 

~^ (^01 ('^mn ^pq ^mp ^nq ~l~ ^mq ^np^ 

— ^Om i^ln Spq — ^Ip ^nq + ^Iq ^np) 
+ ^On {^Im Spq — ^Ip ^mq + ^Iq ^mp) 

— ^Op i^lm ^nq — ^In ^rnq + ^Iq ^ mnj 

+ ^Oq i^lm ^np — ^In ^mp + ^Ip ^mn)] (E.5) 



which is a special case of 



. ' J -J ^2 -J : \ .+ -+- + - 

tr I Cf^ C I Cj^ C nCp C q J = tr (r^Fj rrrtrnrpFg 

4 [z Cklmnpq 

~l~ ^A;( '^mn ^pg ^kl ^mp ^nq ~l~ ^fei ^mq ^np 

^km ^In 5pq ~l~ ^km ^Ip ^nq ^km ^Iq ^np 
~l~ ^kn ^Im ^pq ^kn ^Ip ^mq ~l~ ^kn ^Iq ^mp 

^kp ^Im ^nq ~l~ ^kp ^In ^mq ^kp ^Iq ^mn 
+ hq Sim K.p - hq n ^mp + Skq Sip Smn] ■ (E-6) 

Further traces can easily be obtained from the above results by means of the relation 
[This is a special version of eq. (238) in disguise.] 

-J ; 

Cfc + Cfc = -2 4o U ■ (E.7) 



Appendix F 

In this Appendix we give the explicit relation between the representation of SU (4) 
matrices discussed by Ostlund and Mele [43] (cf. Subsec. 6.2.1) and our parametriza- 
tion (226). The coefficients of the two parametrizations transform into each other 
in 4 groups of 4 coefficients. 



79 



Group 1 reads: 

^ 4 (-^^ ~ ^ -^(-1)0 Ti2 — i T34) , 

^00 = ^{To-i ^(-1)0 + iTi2 + i T34) . 

Xii = ^ (To + i r(„i)o -iTi2 + i , 

X22 = ^ (To + i r(„i)o + iTi2-i Tsi) , 
and the inverse relations are 

To = - (X(_i)(_i) + Xoo + Xu + X22) 

T{~i)o = ^ (X(_i)(_i) + Xoo — Xii — X22) 

T12 = - (X(_i)(_i) — Xoo + Xii — X22) 

T34 = - (X(_i)(_i) — Xoo — Xii + X22) 



Group 2 reads: 



1 



^{-1)0 = ^ (-T13 -iTu-i T23 + T2A) , 

-^O(-l) = - ( Ti3 — i Tu — 2 T23 — T24:) , 

X12 = ^{ Ti3 - i Tu + i T23 + T24) , 

X21 = ^ (-T13 - i Ti4 + 2 T23 - T24) , 



and the inverse relations are 



^13 — ^ (—-^(-1)0 + -^O(-l) + ^12 — X21) 

Ti4 = - ( X(„i)o + Xo(_i) + X12 + X21) 

T23 = ^ ( -^{-1)0 + -^O(-l) ~ ^12 — X21) 

T24 = ^ ( -^{-1)0 " -^O(-l) + ^12 — X21) 
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Group 3 reads: 

^(-1)1 = ^ (-^{-1)3 - i T(^-~i)4 - i To3 + T04) , (F.17) 

^02 = ^ ( ^(-1)3 - ^(-1)4 - ^03 - ^04) , (F-18) 

= ^ (-^{-1)3 + « ^(-1)4 - « - ro4) , (F.19) 

X20 = \{ T(_i)3 + z T(_i)4 - z To3 + T04) , (F.20) 
and the inverse relations are 

^(-1)3 = ^ (—^(-1)1 — ^02 + ^i(-i) + -^20) , (F-21) 

T(_i)4 = - ( — X02 + ^i(-i) — -^20) ; (F.22) 

^03 = ^ ( -^(-1)1 + ^02 + -^i(-i) + -^20) , (F.23) 

^04 = ^ ( ^{-1)1 ~ ^02 — ^i(-i) + -^^20) • (F.24) 

Group 4 reads: 

^(-1)2 = ^ ( + i T(_i)2 + 2 Toi - T02) , (F.25) 

^01 = ^ (-^(-1)1 + ^ ^(-1)2 -iToi- T02) , (F.26) 

^10 = li T(_i)i + z T(_i)2 - z Toi + T02) , (F.27) 

^2(-i) = ^ (-T(_i)i + 2 r(_i)2 + ^ Toi + Toa) , (F.28) 
and the inverse relations are 

^(-1)1 = - ( -^(-1)2 — -^01 + -^10 — -^2(-i)) , (F.29) 

T{~i)2 = ^ (—-^(-1)2 — -^01 — -^10 — -^2(-l)) , (F.30) 

Toi = - (— X(„i)2 + Xoi + Xio — -^2(-i)) , (F-31) 

T02 = - (— X(_i)2 — Xqi + XiQ + X2(_i)) . (F.32) 
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